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ABSTRACT 


In this work, the interface of a non-equilibrium two- 
phase system is viewed as a three-dimensional zone. The usual ap- 
proach, wherein surface fields are postulated instead of being 
derived, is first critically reviewed. A theory assuming conti- 
nuously varying densities in the transition zone, in conformity 
with modern statistical mechanical theories of the equilibrium 
interface, is then constructed. An exact surface balance equation 
is obtained, which is formally identical with that found in the 
usual approach. Complete correspondence with the latter obtains 
in a suitable zeroth order approximation. A hierarchy of balance 
equations pertaining to normal moments of order = 1 is next ob- 
tained. Finally the specific cases of mass, momentum and energy 


are analyzed. 


iV 


aang aie: tdi? view, 7 
ae bad bayer vat ann i 
egy Deas ae pores By) Gak } 
. vey a 


chr as ghee pert oad (| pen 
7 


apts wiak. shea brea theo PS ‘paicenbe: Seah. at 
AT ot Reng) alld ABH Torttncst yotweee) af cote daniel 


ghtepad wate! add AstW sons) ongervos watann irsanace’ fs 
anal sd to tiftareees tit: ms Am itoor ets salvo dtornss sigedive afi i 
+0 deaneet fx aaito: WD si wae RIN OT parritndag ped ; 
eps hay hee sesiteinas ies +4 eon 2h hasge set G hE. Panter: 
beset, eb 


PREFACE 


In the last two decades, many studies have been devoted 
to the non-equilibrium thermodynamics of fluid interfaces. The com- 
mon approach consists in postulating surface densities and current 
densities, and then deriving balance equations relating these to 
the jumps of the bulk fields at the interface. No reference is 
made to the three-dimensional nature of the transition zone in that 
con Cext: 

In erenne hori domain, there exists a well-known 
method, invented by Gibbs, that provides a beautifully simple and 
yet exact interpretation of the surface densities, in terms of the 
diffuse structure of the interface. The object of the present work 
is the extension of this method into the realm of non-equilibrium, 
thereby providing explicit interpretations for the quantities pos- 
tulated in the common approach to the non-equilibrium interface. 

In Chapter I, we give a short review of the Gibbs method 
in equilibrium and then a survey of pertinent theoretical work 
that has been done on the non-equilibrium problem. 

In Chapter II, we first secure some geometric and kine- 
matic tools and then proceed to a somewhat original and systematic 
development of the empirical approach to interfacial dynamics. 
This provides us with a background against which the more general 
theory that follows may be conveniently compared in all details. 


For lack of a better word, and not in a derogatory sense, we use 
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the adjective "empirical" in refering to any theory in which the 
three-dimensional nature of the transition zone is not explicitely 
recognized, and in which surface excess densities are therefore 
postulated instead of being derived. In presenting the empirical 
theory, we strive to isolate a minimal number of basic concepts 
and assumptions, and to derive the rest from them, appealing oc- 
casionally to well established general principles. This enables us 
to pinpoint those results which must be introduced on an ad hoc 
basis. 

In Chapter III, we come to the more interesting part 
wherein the Gibbs scheme is extended into the realm of non-equi- 
librium. We first define appropriate surface normal moments that 
pick up the information necessary for a macroscopic two-dimen- 
Sional description of the microscopic three-dimensional transi- 
tion zone. We then derive an exact surface balance equation that 
is gormaklLy identical with the one used in the empirical approach, 
and explain why this formal identity is remarkable.Next we under- 
take a systematic comparison of the results obtained with those 
used in the empirical approach. It appears that full correspon- 
dence obtains only in a well-defined zeroth order approximation, 
with a wide range of validity. 

In Chapter IV, we derive balance equations for higher 
order moments and give an application to the transversality con- 
dition derived in a (empirical) theory of Bedeaux, Albano and 
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In Chapter V, we analyze the specific cases of mass, mo- 
mentum and energy. While complete correspondence with the empirical 
approach obtains for the first two - at least in the zeroth order 
approximation - discrepancies appear in the case of energy, unless 
some residual terms, of a kind not met in the case of mass and mo- 
mentum, can be neglected. These residual terms arise from the non- 
linearity of some energy related expressions. Interestingly, this 
problem is related to our identification, in Chapter II, of the 
treatment of energy as a weak spot, or less convincing part, in 
the empirical approach. 

Finally a word about the labeling of equations. Equations 
within a given chapter are labeled according to the section in 
which they appear and to the order of appearance in that section. 
For instance, eq.(2.3) means the third equation of Section 2 of 
that chapter. When refering to an equation in a preceding chapter, 
the number of that chapter is added before the other two numbers 
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CHAPTER I 


FLUID INTERFACES IN AND OUT OF EQUILIBRIUM 


1) Review of Gibbs equilibrium interfacial thermodynamics 


Since the object of the present work is the extension of 
the Gibbs approach (Gibbs [1961]) to interfacial systems into the 
realm of non-equilibrium, it is appropriate to start with a brief 
review of that approach. A thorough exposition may be found in the 
book by Defay et al [1966]. It will be sufficient for our present 
illustrative purposes to consider an equilibrium two-phase system 
with a flat interface. The systematic extension of the theory to 
equilibrium curved interfaces has been carried out by Buff [1956], 
whose contributions have been nicely reviewed by Melrose [1968,1970] 
The work of Buff has in turn been extended by Boruvka and Neumann 
[1977]. Curvature effects are also discussed by Sanfeld [1971]. 

For simplicity, we shall first consider a one-component 
system. When surface effects are neglected, all the thermodynamic 
information concerning the system is contained in any one of the 
thermodynamic potentials (U,S,H,F,G) provided these are known as 
functions of their natural variables. For example, knowing the 
Helmholtz free energy F(T,V) we get the equation of state via 


p = -dF and all thermal properties via S = -dF. The liquid-vapor 
OV oT 


transition is also accounted for by F through flat portions of 
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isotherms below the critical temperature and the values of S at 
both ends of these flat portions give the latent heat of vapor- 
ization. Completely missing in this description however are the 
surface effects. But the full picture may be restored by incorpo- 
rating the interfacial area among the extensive variables on which 
the free energy depends. Denoting by U the total energy and in- 
cluding for convenience the mole number N among the independent 
variables, the basic relations are 
dU = TdS - pdV + y dA + LAN 


2 Ue Hs 


dF = -SdT - pdV + ydA + wdN (4S) 
where up is the chemical potential and y is by definition the sur- 


face tension. It follows that 


oF} 
a= ata ee, 


To visualize such a system,one may consider a container whose walls 
may be deformed so as to permit variation of interfacial area A 

at constant volume if so desired. Effects of the gravitational 
field other than that of providing a flat horizontal interface 
(away from edges) are neglected. Far away from the interface, the 
fluid is homogeneous. Extensive properties of bulk fluid regions 
are then given by volume integrals of uniform densities (neglecting 
external fields). Upon crossing the interface or transition zone, 
the values of the densities change from those appropriate to one 
phase to those appropriate to the other phase. Despite macro- 


scopic appearances, these variations are not jump discontinuities 
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but are simply the result of very steep gradients. In simple 
fluids, far from the critical point, the thickness of the transi- 
tion zone is of the order of a few molecular diameters (Lekner 
and Henderson [1978], Goodisman [1979]). The determination of the 
density profile inside this region is one of the main problems con- 
sidered in the equilibrium statistical mechanics of fluid inter- 
faces developed especially in the last two decades (Navascués [1979], 
Croxton [1978]). Experimental investigations based on light ellip- 
ticity have been reported (Beaglehole [1979]), but computer exper- 
iments are the more usual testing ground. 

The diffuse nature of the transition zone implies that 
the partitioning of fluid attributes into contributions from 
both phases and from the interface is not unique. Now, in macro- 
scopic physics, one is not interested in the detailed structure 
of the interface, and Gibbs has introduced a beautifully simple 
and yet rigorous way of keeping track of the interfacial contribu- 
tions, while discarding irrelevant microscopic information. It con- 
sists first in introducing a dividing surface (simply an abstract 
mathematical surface) located somewhere in the transition zone 
and in extrapolating the constant values of the bulk densities 
in either phase up to this surface 2. Inasmuch as the transition 
zone, despite its extreme thinness, is diffuse, the computation, 
with such extrapolated densities, of the total amount of an attri- 
bute in a region V encompassing a patch from the interface will 


usually be in error. The second aspect of the Gibbs scheme consists 
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in correcting these errors by compensating surface excess quantities, 
which are the surface integrals of superficial densities defined on 
the dividing surface %. This method is very attractive and widely 
used because it is in tune with macroscopic intuition and yet ex- 
tremely general and exact. This last feature is not always appre- 
ciated as the scheme is now and then refered to as the "Gibbs approx- 
imation" in the literature. 

The gist of the approach can be clearly illustrated 
through Gibbs' definition of adsorption. Consider a two-phase fluid 
mixture with bulk mole number densities pees a » where the super- 
script a = +,- identifies the phase. Le us denote by N’(V) the exact 
mole number of species < contained in a cubical region V with two 


faces parallel to 2. The dividing surface splits this region V into 


two subsets V" and U™ and N’(V) is written as 


N*(V) = | n* dV + | n. dV + | PedA false, 
AG ee 4 A 
V V Gy) 
Saye ie es ee 
a A A 
where X(V) = Zn V_ is that part of © contained in V and r, is by 
definition the adsorption of component <. It is an intensive sur- 
face variable. Denoting by c the coordinate on an axis perpendi- 
CGuiate tOure (ceo Oro »)eandepy n the exact mole number densities, 


which vary continuously across the transition zone, it follows 


immediately that 
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The diffuse nature of the transition zone confers a dearee of arbi- 
trariness to the choice of the dividing surface, and the values of 
the surface densities usually depend on this choice. This dependence 
on the location of © is easily understood in the case of the adsorp- 
tions defined above. Let us assume for example that ne < ni » and 
that n(z) is monotonous across the transition zone. If the dividing 
surface is located at a value of Zt where n’ has reached the value Le 
it 1SeCléeanethat r; will be negative because the sum of the integrals 
of the bulk densities then overestimates N(V). On the contrary, put- 
ting at a value of zg where n has reached the value ni leads to an 
underestimation of the mole number and thus to a positive value of rie 
This case provides a dramatic illustration of the possible sensitive- 
ness of surface density values on the location of the dividing sur- 
face: by moving © within the extremely narrow margin delimited by the 
transition zone, the adsorption re can even be made to change sign : 
There are however noteworthy exceptions to this type of behavior. 

Let a denote the exact density of some attribute, with corresponding 
bulk densities a’ and a. It is clear that the associated surface ex- 
cess a> will be independent of the dividing surface position if AS 
a. One case in point is the surface tension of a flat interface 
(Green £1969]). Another example of almost complete invariance is that 
of a dilute solution of a strong surfactant (Lucassen-Reynders [1976], 
p.282). This invariance explains why surface chemists can refer loose- 


ly to the adsorption as "number of molecules per unit area of ad- 


sorbed film" without running into inconsistencies. 


* a , 
; 


OV 


° S. at baptert a Be iusche caw oe 
te  Frie Y | 1 Renin se Sama gz tut “lsvuiie Ban Mab 
een it sts ‘hoksrenesy aul peo eupite) Oya at \ih 
-iouTay wilt Darou ‘ang seh vin) pv Ahoy, 9 sia, behagny ef so 
Coca eal oe ae att Seva Oy set td Gain 7 SBR? «pata-2h 
oa Sadana a thir epider’ 4237900 fight aut stag ud a, # 

he. oF elise} en sae sis eadsedd abi vt 630% 3.90 stiri & Sa, ae 


> 


my Va seth eee i bie Suit fine 43 es pS Faia ooh te nut ane eet 
soem ach 09a and ty oa Saint Vie Ree ateyeas Aebs BAF 
wile seibiad” ate 70 ry okewor Siiong balay rary at he th aren 
nits Det thattah sip» nosIAN Ulowey tes att -niitwid anion amt" 
> arenidilt: 6a sbent’ ed nove As, CORRE DH BF 0 Ot aaty recon 

eh upist 40 eay2 site ot ant REOKE LONG TY sult wen anne 
prides hi0s AS iw: yaa 9 1 pie 7 Ww ; ani) tags sttopitdnes “se tod 
ae Soa Ie bade dee, Ey a i Sqolorei ol ce Snes 2dberenet sad 
oh nothte09 naira fe 09 sit! Yo teghnugahints 4d [lew ‘ome 

| erat iont yor Lf ‘ aan 6 #94; tye sq: at Yahog ni geaa Bee a 7 
font af soem: sri saowhs ohne TSO BT) ¢ ) 
ewer ; rab Anutsaewe ony ite ci%e neuen suite we | | 
raeebe lassie arn 217) tts adit LABY ute | 


| Osea” aD ms 
a pitennv', Uae a | At 
> 


Even in cases where the surface densities depend sensibly 
on the choice of the dividing surface, the theory can be fully devel- 
oped leaving that choice open. The choice of a particular dividing 
Surface becomes then largely a matter of convenience and varies with 
the context. Usually specific choices are formulated in terms of 
adsorptions, for instance by requiring the adsorption of a certain 
component to vanish (Defay et al [1966], p.27). 

To further illustrate the working of the equilibrium theory 
we shall briefly consider the surface Helmholtz free energy in a one- 
component system. Once the values of T,V,A,N are given, assuming 
that some convention (which need not be specified) has been made on 
the choice of the dividing surface £, the total free energy is writ- 
ten as 

Elis Veteiiim=e ro (lave, Neen ns (1 VaGNmaee re Gl VRAGN): ORX1.5) 
where the values of V',N’,V~.N7 depend on the choice of &. The free 
energies F’,F refer to the two bulk phases and F° is the surface 
excess. The pressure in the system is known and so are the mole num- 
ber densities n’ and n~ in the bulk regions. Now let us, at fixed in- 
terfacial area, change the mole numbers in the bulk phases by adding 
isothermally amounts AN” and AN and changing the volumes so as to 
accomodate this new material without changing any of the intensive 
variables (whose values, according to the phase rule, are all deter- 
mined by T along the liquid-vapor coexistence line). This may be 


achieved by taking AV - Lant and eA ee ea . The new overall 
n 


n 


volume and mole numbers are denoted by V' and N'. Clearly the varia- 
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tion F(T,V',A,N') - F(T,V,A,N) will be fully accounted for by the 


variations of the first two terms on the right-hand side of (1.5). 
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where, since ie n , the above 2x2 matrix is invertible. It follows 


We have 


that V and N can be varied independently without affecting the value 
of the last term in eq.(1.5). Whence F° is actually independent of 
V and N and depends only on A and T. Furthermore, it follows by homo- 
geneity that at fixed T, F° is proportional to A, that is F°(T,A) = 
Af>(T) . We may thus finally write 

FCI VN Aah (tov Neb a(TeVeeNGagAt Gili (126) 
We emphasize that the explicit form of this function f>(T) depends 
on the choice of the dividing surface. As we shall see below, its 
physical meaning is most transparent when the dividing surface is 


chosen so that the adsorption, which in the present case is given by 
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vanishes. We note incidentally that the left-hand side of eq.(1.5) 
makes no reference whatsoever to a dividing surface, which im- 
plies that the surface tension (see eq.(1.2)) is independent of its 
location. In principle we could work with this thermodynamic poten- 
tial F only. But the decomposition effected in eq.(1.6) is more ef- 
fective in assessing surface behavior. The gunctions F* and F~ are 


also independent of the choice of £, but V", V, N* and N7 are not, 
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which is why the values of the first two terms on the right-hand 
side of eq.(1.6) do depend on the location of ©£. 

A number of thermodynamic relations may now be easily ob- 
tained. We first decompose the energy and the entropy in the same 
manner as the free energy: 
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SCTRVER INS OTS V aN (er eSaGns VeeNe) ERS CEA) 
so that 
pat eearse (1.9) 
Differentiating eq.(1.6) we get 
Srsidue- ody, hes dN =e 9S dy = pay 4 ed” 
eed) epee die de 
When combined with eqs.(1.8) and (1.9) this gives 
dF> = - S°dT + y dA + pdN° (1.110) 
du° = TdS> + ydA +udN> . eit) 


One must be careful not to misinterpret these relations: T,A and 
N° are not independent variables. 
Integrating eq.(1.11) at constant temperature yields 
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from which, using eq.(1.9), 

es vy A + ne (epae) Aue Cy tee) 
Upon differentiating the last relation and comparing with eq.(1.10), 
we obtain a surface analogue of the Gibbs-Duhem relation, known 


as the Gibbs equation: 


dvaee=as: din tednis (1.14) 
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When the T = O convention is adopted for the choice of the dividing 
surface, eq.(1.13) reduces to 

feee yoo (1.15) 
This is in agreement with the fact that y depends only on tempera- 
ture in a one component system. We note that the identification 
between y and surface free energy per unit area obtains only when 
the T = 0 convention is adopted. Discussions of the temperature 
dependence of y may be found in Defay et al [1966], Gugenheim [1965], 
Buff [1968] and Bikerman [1970]. Equation (1.14) also provides an 
interesting relation between y and the surface entropy per unit 
area in the [ = 0 convention: 
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Naturally S° is related to thermal effects. Let us consider a change 
of the area A at fixed T,V,N. In principle, a variation of A at fixed 
T and N implies a certain redistribution of matter from bulk to in- 
terface, and may thus entail some variation of V. But, unless a very 
large amount of area is created, such a change will obviously be so 
small as to be completely negligible. The heat absorbed by the sys- 
tem in the course of such a reversible isothermal expansion is given 
by 

GOM=t ads 


N 
denotes the total entropy at fixed N. We thus have 
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Together with eq.(1.16) this implies that, with the T = 0 convention, 
0) 287 (1.18) 
dA T 
In the general case, this is replaced by 
| i | =T(s?+ PGE) (1.19) 
i dT 


All the above relations in which the [f = O convention is not 
adopted may be generalized in straightforward fashion to multi- 
component systems simply by replacing terms involving T and yu by 
sums of terms involving the individual adsorptions Pr; and chemical 
potentials u, (Kirkwood and Oppenheim [1961]). 

The equilibrium statistical mechanics of the fluid inter- 
face has undergone significant developments in the last two deca- 
des (Navascues [1979], Croxton [1978]). It makes extensive use 
of reduced distribution functions and is based on the idea of 
continuously varying exact densities in the transition zone, 
whose meaning is based on probability. The theory of Chapter III 
rests in part on the assumption that such meaningful densities 


also exist out of equilibrium. 
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2) Studies of non-equilibrium interfacial systems 


In this section, we present a brief survey of studies 
devoted to the non-equilibrium fluid interface in an attempt to put 
the present work into proper perspective. In view of the over- 
whelming abundance of related literature, we shall restrict our- 
selves to a narrow path that will lead us to the more recent stu- 
dies to which the present work is connected. For reviews, the reader 
is refered to Levich and Krylov [1969], Lucassen-Reynders and Lucas- 
sen [1977], Defay,Prigogine and Sanfeld [1977],0strach [1979]. An ex- 
tensive bibliography covering older work, particularly in connection 
with free thin liquid films, may be found in Mysels [1959]. 

The most well-known phenomena in relation with interfacial 
motion are probably the Gibbs and Marangoni effects (Kitchener [1964]). 
The Gibbs effect is in reality a static type of elasticity operative 
in thin liquid films (Krotov and Rusanov [1972]). When a thin film 
of a surfactant solution is stretched, the interfacial area on both 
of its faces increases. As a result, additional surface-active solu- 
te migrates from the internal solution which is thus depleted with 
respect to the solute. In the new equilibrium that is established, 
it follows that the surface tension is higher. Compression of the 
film naturally produces the opposite effect. It is obvious that this 
effect must depend on the film thickness. A thermodynamic theory has 
been given by van den Tempel,Lucassen and Lucassen-Reynders [1965]. 

The Marangoni effect (Scriven and Sternling [1960]) is re- 


lated to the Gibbs effect but does not hinge upon the thinness of 
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bulk substrate. When the interface of a solution of a surface-active 
agent is rapidly expanded, there results an impoverishment in solute 
near the interface, which is accompanied by an increase in surface 
tension. The reestablishment of equilibrium is delayed by the dif- 
fusion of solute in the bulk. This explains why a transient higher 
surface tension may exist without the requirement of a thin bulk 


phase, as opposed to the Gibbs effect. The best known early theo- 


retical investigation in that field is that of Ward and Tordai[1946]. 


Further significant developments were made by Baret [1968] whose work 
contained a mistake later corrected by Petrov and Miller [1977]. A 
modern review of dynamic surface tension has been given by Defay 

and Petré [1971]. 

Another important concept in relation with the fluid in- 
terface is the so-called surface viscosity. This concept was intro- 
duced apparently for the first time by Boussinesq [1913] and later 
formalized by Ericksen [1952], Oldroyd [1955] and Scriven [1960]. 

The theory accomodates shear as well as dilatational viscosity 
(Slattery [1978]). Non newtonian behavior has also been considered 
(Addison and Schechter [1978]). Reviews have been given by Joly 
[1964] and Goodrich [1973]. Another phenomenon that has been stu- 
died is surface diffusion (Brenner and Leat [1978] and refs therein). 

In 1959, Sternling and Scriven made the first attack on the 
problem of interfacial instability and turbulence. Interest in this 
problem has been recently revived (Hennenberg et al [1980] and 


refs. therein). 
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Certainly a landmark in the subject of interfacial dyna- 
mics was the paper by Scriven [1960]. It was followed by a series 
of investigations by various authors leading to recent systematic 
applications of non-equilibrium thermodynamics . His work made no 
explicit reference to the Gibbs approach and was in the spirit of 
the empirical approach discussed in Chapter II of this thesis. In 
a first step, Scriven took an intrinsic view of the interface with- 
out reference to the adjacent bulk phases and wrote down equations 
of motion within this two-dimensional world. In a second step, he 
considered the connection with the bulk phases. A pedagogical and 
detailed exposition of the scheme invented by Scriven has been given 
by Aris [1962]. 

Scriven derivations however entailed some conceptual 
difficulties which were pointed out and corrected in a generaliza- 
tion taking proper account of mass transfer across the interface by 
Slattery [1964]. 

Next came a paper by Ghez [1965] containing two innova- 
tions. First, the generic surface balance equation (at the empiri- 
cal level) was derived; secondly, it was the first time (to our 
knowledge) that the standard methods of linear irreversible thermo- 
dynamics were brought to bear on the problem (though it was only 
later that they were fully exploited). Building upon the work of 
Ghez, Georgescu [1969] considered momentum balance at an interface 
with internal angular momentum. Ghez also made specific applications 


of his theory (Ghez [1970]). 
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Waldam [1967] also applied the methods of irreversible 
thermodynamics, though without considering surface energy, to the 
boundary conditions at the interface of two immiscible fluids. His 
work was generalized in a paper by Bedeaux, Albano and Mazur [1976], 
who included surface energy thus bringing in surface tension. These 
authors set up a formalism (hereafter refered to as the BAM theory) 
that has since then been used in several other investigations. 
Their work has been extended first by Kovak [1977] and by Vodak 
[1978 b] to a multicomponent system with what they call singular 
mass densities at the interface. The entropy production at the in- 
terface in such systems has a rather involved structure and there 
results a rich supply of phenomenological coefficients. The theory 
was further generalized by Wolf and Albano [1979] to accommodate 
electromagnetic fields. 

We should also mention another work of more limited scope 
by Popielawski [1970] and a careful treatment of diffusion related 
matters by Kehlen and Baranowski [1977]. 

Other approaches, in the spirit of rational thermodynamics, 
have also been made (Moeckel [1975], Murdoch [1976]). 

A common feature of the above theories is that,in their 
context, surface densities are postulated instead of being derived 
and that they at best only implicitely acknowledge the three-dimen- 
sional nature of the interface. Interfacial balance equations were 
investigated by Slattery [1967] from a more fundamental point of 


view according to which the transition zone is three-dimensional. 
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This is in agreement with modern equilibrium statistical mechanics 
theories of surface tension and related matters (Navascues [1979]). 
Slattery presented his work as a generalization of the work of Buff 
[1956] on equilibrium two-phase systems (further developments in 
this equilibrium theory were made by Boruvka and Neumann [1977]). 
In the present study, specifically from Chapter III onward, we 
start from the same basic premises that were neatly formulated in 
this paper of Slattery and carry the theory considerably further 

so that all the aspects of the empirical approach may be interpreted 
in terms of appropriate surface normal moments. Some of our results 
have been reported in a paper which will appear in Physica A. The 
same problem has been approached differently by Ronis, Bedeaux 

and Oppenheim [1978] and by Albano, Bedeaux and Vlieger [1979]. An 
approximate derivation of the surface balance equation was also 
given by Deemer and Slattery [1978]. 

In the next chapter, we shall look at interfacial balance 
equations from the empirical point of view, that is without refering 
to the diffuse nature of the transition zone. This will provide us 
with a background against which the subsequent theory may be com- 
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CHAPTER II 


THEORIES POSTULATING SURFACE DENSITIES 


1) Geometric preliminaries 


In this section we set up the notations used throughout 
this work, discuss some key geometrical concepts and derive some 
kinematic results needed in the subsequent discussion. The first 
part is concerned with the calculus and geometry of surfaces and 
the second part is concerned with their motions. We shall view a 
Surface as a two-dimensional submanifold of the three-dimensional 
Space E3, with or without boundary. Surfaces with self-intersec- 
tions are thus excluded. 

The tangent plane to a surface © at a point p of %1s 
a two-dimensional vector space denoted by Ta(2). To indicate that 
a vector is tangent to x, we will put a bar over it. Any vector 


eo Es may be decomposed into its normal and tangential components 


(X) = 0) ’ (Ags) 


where Pe is the projector on ie and n is the unit normal vector 
to © at p. The directional derivative at Dp, with respect to 
Ue Ta(2), of a function f (resp. a vector field X) defined on 


is defined by 
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where — + C(E) is any curve on & such that ¢(0) 


tl 
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I Rice) is a parametrization of = around p « Z, and 


U = utdR + u29R , then 
Paes 0E2 


1 
f= 
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Df (P) 


and likewise 


DX (B) = uta__ (FoR) (B) + u2a_ (RR) (B) 


where the symbol "o " stands for the composition of functions. Note 
that U need not be a unit vector and that both Df (p) and D{X(p) are 
linear in U, at fixed p. Thus the vector Gradf (p) « T(z) is defined 
by 

Do f (p) = U-Gradf (p) for every U « Ta(). flee) 
The covariant derivative at p, with respect to Ue Ta), of a 
tangent vector field X defined on £, ¥5X(P) » is defined as the pro- 


jection on the tangent plane of its directional derivative: 


gate (1.4) 


> 


V,X (Pp) = is 
At fixed p, the map U X(p) is clearly a linear operator on 
Ta(2)s it is denoted by (VX)> and called the covariant derivative 
of X at p. Hence by definition, 

[CV X]+(u) = V_X(p) for every U e T(z). (1.5) 
Given a linear operator A : T(z) > Te), its trace will be denoted 
by Tr(A). (The bar emphasizes that we are dealing with a two-dimen- 


sional vector space.) Given a tangent vector field X on £, the diver- 
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gence of X at D e © is defined as the trace of its covariant deriva- 
tive at D: 

DivX (p) = Trl VXI. (1.6) 
Let A be a field of linear operators on £, which to each D asso- 
ciates a linear operator Ae that maps TAC) into itself. The cova- 
riant derivative of A at D, with respect to Ue Ta), is defined by 


(V-A) = Pal Dane): (restricted to Ta(2)) : (ez) 


ae 
p 
At fixed p and U this is, just like at a linear operator on Ta(z) 


and it is easily seen that (V5A)> s (V5A" )s , (where AX denotes the 
adjoint of A ) and that VoCA(X) J = (VoA)(X) - A(V_X) , where X is 
any vector field on £. Finally the divergence of A at D is a vector 


DivA (p) « Ts(Z) defined by 
A = a * _ = aS =, 
DivA (p) = CV, Adg (Ua) ~ Cote (ey = (1.8) 


where (U;,U,) is any orthonormal basis of Ta(2). If X is a tangent 
vector field on x, the following identity holds: 

DivlA(X)] = X-DivA + Tr(A VX) (1.9) 
where VX is the covariant derivative of X (see eq.(1.5)) 

On each tangent plane Ta) of a surface © is defined a 
linear operator s Tale) - Tal) called the shape operator of © 
at p (O'Neill [1970]), which tells how the unit normal n to © chan- 
ges in the neighborhood of D. Speciiicallys 1h Ure T(z), we have 


on +> ),> 


R(B) + Deh (p) = 40 


ae 


since nen is constant (=1). It follows that Dn(p) E Taz). The 


shape operator at p, Si , is the linear operator on Tale) defined by 
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Ss (u) 25 D5n (p) for every U « Ta). aD) 
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(The operator 0 Fs is usually called the Weingarten map, whereas 
the twice covariant tensor associated to % is called the second 
fundamental form.) It may be shown (Hicks [1965]) that the shape oper- 
ator is self-adjoint, that is 
(u)*V for every U,V « T(z) : (Teste) 
Given a unit vector @ « Tal) the plane through Ds Spanned by n(p) 
and G, cuts ¥ along a curve T,called a normal section at p. The nor- 
mal curvature of ¥ at D in the u direction, Kp(G) » is by definition 
equal to the curvature of the normal section TS at p, counted as 
positive when ly bends toward n(p). It may be shown that (O'Neill 
[1970}) 

Kp (G) = G-S3(u) : el) 
Since 3 is self-adjoint, it has two real eigenvalues k,,k,, called 
the principal curvatures of x at p. From eq.(1.12) these are the 
extreme values of the normal curvature at p. The mean and Gaussian 
curvatures of © are defined respectively by 

H = 4(k,+k,) = 41rS and K = k,k, = detS . Gieey 
The sign of K(p) tells much about the shape of © in the neighborhood 
Gf p (O'Ned 11) £19709.p.204). 

life dee Ta(z) and X is a tangent vector field on ¢, 

neDoX = =X-Pen = X-S(U), a result which combined with eqe(1.4) im— 
plies the Gauss equation 

DX = ¥,X+ 0-S(X)n. (1.14) 
An immediate consequence of this equation is that if A is a cons- 


tant vector and A is the tangent vector field defined on 5 by 
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ws 2 hats. (1.15) 

Finally we consider a moving surface. This will be viewed 

as a one-parameter family of surfaces Dye one for each time t, such 

that there exists a reference configuration £° and a differentiable 

function R(tsp); D « ©° ,with the property that the map Ree Die 

pe R, (6) = R(t,p) is a diffeomorphism. The reference configuration 
>° need not coincide with any of the surfaces of the family. More- 

over given the one-parameter family re and the reference configura- 

tion 5°, the parametrization R(t,p) is not unique. We will return to 


this point shortly. The normal speed v- of the moving surface at 


ie Ge) e Eisedetined iby 


dP (t) = viA(p) (1.16) 
is 
wan 5 - wae = ace 
where P is the trajectory defined by P(t) = p, P(t') « Da and 
se ae : : 
nh edP(t) = 6 for every t'. If a function f is defined on the 


dt 
moving surface (that is a function i is given on each r,) the nor- 


mal time derivative of f, Celie. is defined by 
dt 


d'f = d £(P) ea) 


where the trajectory P is as above. Given any trajectory q(t) such 


that g(t) « Z, for each t, it may be shown (Treusdell [1960]) that 


t 
neq =v (1.18) 
and that 


“Grad f = Gy (1.19) 
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If a tangent vector field X is defined on the moving surface and 
Xe is the vector-valued function of time Xa(t) = (XM) eeand eit 


(@,,€,,¢3;) is a fixed orthonormal frame, then 


Wy Ges es 


where use has been made of eq.(1.15). Together with eq.(1.19) this 


q, -Grad (@, +X) SS OWES Jos NOTES = 6; [VEX + esta 


implies 


Ap ic VERE Ges One (1.20) 
dt 9 dt as 


We shall finally derive expressions for the rate of 
change of the unit normal and local area of a moving surface in 
terms of its normal velocity field. Given two surfaces © and 5', 


and a smooth map M : Y + r' between them, the differential of M at 


Dp « Z, denoted by The is the linear map between Ta?) and Tp. (2") 
= =e ° 
(p' = M(p)) defined by 
Hie (Z) = OLM(p) for every U< Ts(z) . (1.21) 


Given p « ©, the jacobian of M at p, U(p) , is defined by 

J(p) = n+ (MCU, XM a(t) 
where (U,;,U2) is any orthonormal basis of T(t) such that 
n-(U,xU,) = 1 and n' is the unit normal to x' at p' = M(p). The 


intuitive meaning of the jacobian is revealed by the change of 


variable formula : 


> 
| f dA == | foM |J| dA Cip2s) 
¥ ' 
where dA and dA' are the area elements on © and y' and f is any 
integrable function on r'. 


Coming back to a moving surface %,, we assume that its 


c 


a 2 : : 1 © | ; os 7 | 
niece Santee Eg Ml ae 
intomiferyiren tigi tee ae 

a), ae 2 vith, Piatiralteel es ‘Gar Wikhon tia) st” v0 gti. 

. Fibs Sona i saved .bfa'* grtantoy | met edt he emt 
Sr) Hpien yt sea eeh fe t=) TF aa kagizoe > 


Wels us pie. Yet hee ad bp her Fadl “ Nat teSGhan. 2 at 
r a iwt- §} 
Ree . 


WarBanrtab 2! ee ci br i ae Wht aoset e%2 32 iF neta. 
: mits ey Oe Me ‘ . ft 
teat. Hoye fades *, feaal Qmrnmhontere vise 2¥ ally Aly oie 
att Gi 18 te,’ ‘f AN Lain ve finu gd ate iN baie pis (ds sa: 
fu Sphono ody eh palhavnr at metiddey 42 “o gateebont avy Enasel 


aed 


7 6) Want ei aiehtly 
co Shah! et ce ae | 6 
’ ra WS ™ a ; 
; a uth, at > bon Sabin ho eiptnele we of) ars “sale ait 


lng: tess baet ol 


Ge 


motion may be described with the help of a reference configuration 


Oe Z, and of a smooth function R(t t,he 2” eausuch that. alt 
0 


eS 
fixed t, the map R, phe aes Ben R, (6) = R(t, p) is a diffeomor- 
phism between 5° and re: This brings up a subtle point with inci- 


dence on the theory that follows. If the moving surface a were a 


model of say a membrane, there would be a natural choice of the 
function R(t,p) as a "material" parametrization, in the sense that 
R(t,p) would be the position at time t of the membrane "particle" 
which occupied the position D in the reference configuration 

yes Ue at time t,. But when ty is of an abstract nature, with no 
"material" existence, as is the case with the dividing surfaces 


considered below, there is no such natural choice. In fact it is 


clear that given a moving surface ©, and a particular parametri- 


c 
zation R(t). we can find infinitely many other different para- 


metrizations simply by superposing internal "motions" within the 


moving surface. However, given the one-parameter family x,, all 


t 
these different parametrizations have in common the normal com- 


ponent of their velocity field for, according to eq.(1.18), we 


have 


where v- is defined by eq.(1.16). Given some reference configura- 


Clones this leads us to define the normal parametrization R(t, p) 


ty 
based on De as that parametrization for which each point of the 
0 


moving surface has at each time its velocity perpendicular to the 


surface on which it lies at that time. This is tantamount to the 


7 : ea _ = a 
216 wr re pied acbuity pw: pa : irae» dia ahi 
ve yk 2-01 ieee se TCAD ‘7 


p 


aint 2T sant a te ya oR yr 2 fh qed 3 ti 


tag aketa anton efadla « qt etutedatFh y yt baw *r ods a 
A wt", fi Sosaribe vrenvon sh 1}. 2pofbo? geil? yy ioast Say na aonb 
wat Ie Borodo -Tewion 6 4c iow. 27s? .s le YaNSe! & “se wn | abeom 

' a” -_ 


4 


yey ezise one Al. fer fastsigiensg. "TASTE Bye? sae )F net gat 


1 adie ‘si 
P3f9 " Chi May S07 re wird) je OTe Teoy an SO ore (4,2 


a arya 
‘ . b ; ; - 
al SORROTIAGD She HG aes nd oo there eng bo Dgeeepee tte 

3 | 

é 
_ 
. - ‘ « - - 4 
trvirttew, eaurey PSeV sone WR ey .s nore Si ot OT Ph, a Z 
¥ 

a ‘ j { 7 A~< 4 ’ e ul 
Setarie HhPpevre. ant oiew Jets fy ¢ 2 er i - é»y “J 3ST%oy am 


ot 4) Seo of). oatoly. Pete wi tioue! on) stadt wwfed “veh znep 


' . ; - . D ra 
< ? ¢hanG% é{ueltaa ¢ bre. c Baste aitvem & tavea C847 Shee 


=-S780 QNs1eT TES YeAt ee Ylovtn! TT aG> (BWP iS a neeD 
md bet iw ‘sido’ Snyett | pry Pty roqity. Uh. yh? > <r AvtToM. 
, mis 

fin txt ins? tare or? po ne ey" “tl =~. 405 ue CRP wa 
| 

wes, Tena 982 vega at tvet réhe ; at eesti sgene 


oy (EL1i.pa Ot palefonse Solr be) * 7 20)s~ “hee Shee 


- ) 2 then 
> os <6 
rina M09 aanene ; Srg2 nayha. Loh= ap ton 
= a —_ _ 


VA nities intan phan oon SatVi of ia alt Ben 
| > a 
mn inte ee abs | eee Nae 
_ Runs 2 Fuge a ay eye non apt ten a +. 


= : a v 
. ‘ = 2 on + . . ve 7 
acy 6G TRECs eo | a obe C1 unr Wi 
7 = © -¥ 7 . 7 » r 7 _ , 


and that aR* (tsp) 
t 


requirements that R(t, sp) = D for every D € Ly 
0 


Q 


be perpendicular to z, at Ri(t,p) for every t and p « De (Natu- 
rally in general R (tp) z p when t z t,.) We will now study the 
rate of change of the unit normal and local area in the normal para- 
methization. For Simplicity we put ty = 0. Let ue Tp (20) and define 
a "normally convected" tangent vector U(t) by 

U(t) = (ekg (@) (1.23) 
where (Re dep is the differential of R at the point p (see eq. 
(1.21)). Given a coordinate system (£!,¢2) on a neighborhood of p 


in Zo? with u = uraR (0,6) + u29R (0,p) , we have 


de Oe? 
U(t) = utoR (t,p) + u2aR (t,p) 
ia gE? 
SO that 
dU = uta aR + uta ak 
dt 9619 9679 
from which 
dU = D(veh) = weGradven - v'S(u) (Gie24)) 
ae || eae 


where use has been made of eq.(1.10). This is the basic relation 
from which local rates of deformation in the normal parametrization 
may be extracted. Let(U, ,U,) be an orthonormal basis at p with 


NM = U,xU, and U;,(t) and U,(t) be defined by "normal convection" of 


i 
t 


at R(t,p), we have 


Ue. late) Ws aks (Ge) ieee (62S) 5 the Gat) 1S) the jacoblansot R 
at ieee gnand na(t) is the unit normal to De 
na(t) = J) (tb) Ui(t)x Up (t) 


from which, using eq.(1.24) and J, (05p) = 1, we get 
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d nm = he n = ve S(U, )xU» t U,xS(Up )J 
dt 7 |t=0 at | t=0 
+(U, Grad ‘igh nxUs + Us «Grad ve Uixn 
= = |ddu # v- 2H \n.= Gradv- 
OG | tev 
Since d n is tangent to 59, it follows that 
dt t=0 

En = -2Hv~ (dyes) 

a taeD 
and dn = =Gradv . (1.26) 

dt 


These two formulas are remarkably simple. The first one gives the 
instantaneous rate of change of area of a small piece of surface 
whose points are propagating according to the normal parametrization. 
This follows from the geometrical interpretation of the jacobian 
derived from eq.(1.22). The second one relates the rate of change 

of the unit normal along a "normal" trajectory to the normal velo- 
city gradient in a most direct and intuitively satisfying way. 


Finally if r(t) is any trajectory such that r(t) « Z, for 


-)> 


every t, and m(t) is the unit normal to zr, at r(t), it follows 
from eqs.(1.19) and (1.26) that 
di = # Grady = 56) ) where r, = P(r) . tS era 


Hereafter the symbol "Grad" will be used exclusively for functions 
defined on © and the symbol "Div" only for tangent vector fields 


on =. in the sense of eqs.(1.3) and (1.6) 
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2) A minimal empirical theory of the generic surface balance 
equation . 


This section will serve as a reference frame for subse- 
quent discussions. Its purpose is to establish a general surface 
balance equation for a two-phase fluid system at a completely empir- 
ical level, with the minimum number of assumptions. The essential 
results are summarized in eqs.(2.19)-(2.24). The eventual origi- 
nality of this derivation lies in the economy of assumptions and 
in the obtainment of a functional form of the line flux (eq.(2.18)) 
which, to our knowledge , has not been reported before. 

It is assumed that the interface can be modeled by a 
region of zero thickness, a true surface , on which surface den- 
sities and current densities are postulated. The theory is termed 
"minimal" for two reasons. First because we do not adopt any spe- 
cific interpretation for the surface quantities in terms of the 
diffuse structure of a real interface.(This will come later.) 
Secondly because we attribute a priori no tangential velocity to 


the surface ©, representing the interface. But of course, just as 


L 
with any other surface, the moment we assume that it exists (be 

it only conceptually) and that it moves, it follows from the 
considerations of the preceding section that at each time it has, 
at each of its points, a uniquely defined normal speed v- . On the 
other hand, to accept a priori a notion of tangential velocity 


(with presumed physical meaning) for the points of a would amount 


to attributing more or less permanent labels to them. Although ob- 
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viously nothing prevents such labeling in a purely formal or mathe- 
matical sense, in order to have physical meaning it would have to 
be based on some physical picture such as te being material or being 
a membrane. But since our aim in this section is precisely to iso- 
late or define a general empirical framework free of any such special 
assumptions, we must refrain from using the concept of a physical 
tangential velocity of De as a primitive concept. So it will not 
enter into our formulation of the basic assumptions below. 

Because of the small number of assumptions, the results 
derived in this section will be quite general. They will serve as a 
kind of frame of reference for our discussion of existing phenomeno- 
logical theories in the next sections,as well as for the develop- 
ment of the theory of Chapter III, where a specific interpretation 
of the surface quantities will be adopted and investigated in con- 
siderable detail. 

Consider a non-equilibrium two-phase system. The inter- 
face (at time t) is represented by a surface Dy coinciding with the 
phase boundary which will be assumed to be macroscopically sharply 


defined. The surface ©, divides the material into two regions My 


t 
and My pobothe containing Dt for convenience, that is My n M, = 
(the symbol "n" stands for the intersection of two point-sets, 


Re 


that is the set of all points that lie in both sets) and the unit 
normal n to De pointing towards My . According to the continuum 
or "smoothed" description of matter, the fluid is endowed with 


various extensive attributes,the amount, flow and internal produc- 


tion of which may be computed using appropriate densities, in the 
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bulk regions. If V is a region fixed in space not inter- 
secting Z,, the total amount of attribute A inV,A(V) , its flow 
©(3V) through the boundary aV and its internal production rate P(V) 


are given by 


Atti = Jag (2.1) 
V 

o(3V) = | U°N dA (22) 
3V 

p(y) = fo ay (2.3) 
V 


where a, Jando are the bulk density, current density and internal 
production density respectively. These three quantities will be 
hereafter collectively refered to as "densities". At the interface 
the densities associated with some attributes may have disconti- 


nuities in the form of jumps. We will denote by na (resp. a) the 


(resp. M, ) and equal to a in My (resp. M, ) 


function defined on Mt A 


G 
deprived of ry and to the limiting value of a from the + (resp. 
the — ) side of Dee at each point of rye When the density a has a 
jump at a point p of Bes 
trary and in any case is physically irrelevant. We shall adopt the 
convention that a(p) = a'(p) at every p« Z,. The same notations 
and convention will be adopted for J and o. Following the notations 
of Vodaék [1978], the jumps at the interface will be denoted by 


double square brackets: 


fan) = a (pe) =a (p) 
Tol (b) = of(B) - o7 (Pp) (2.4) 
adnate (oer) 


the choice of the value a(p) is quite arbi- 
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for every p eX Next we postulate that, when a {ixed region V 


t 
intersects the interface, the total amount of attribute A in V, de- 
noted by A*(V), is expressible as 

A*(V) = | a dV + | a> dA (2.5) 
V 2 (V) 


where a is the same bulk A-density as above, a> is a surface excess 


density defined on 4. and z,(V) = Vn Z, is the patch of £, con- 


t 
tained in the fixed region V at time t. Equation (2.5) may be con- 
Sidered as the general expression for the total amount of A, valid 
whether V intersects © or not, if the surface integral is set by 


definition equal to zero when V n xX, is empty. Likewise the rate 


E 
of internal production is assumed to be expressible as 


aan Ca ee | odV+ | o> dA. (2.6) 
V D4(V) 


Finally we postulate a surface current density j defined on De and 


such that when a 4ixed surface S, with unit normal vector field N, 


meets Zr, at right angee all along the curve Snt, then the net rate 
Of transter of A through :S 1s 
o*(S) = | JeN dA + | j-N dd (267) 
& Snzy 


where J is the bulk current density and ddA is the arc-length ele- 
ment on the curve Snz,. If the surface S meets re but not at 
right angle, then we write 


oS) = | 3 ee [G9 PAR & (2.8) 


S Sod, 
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where at any D € Snz., Vv is the normalized projection of Noon Ty) 
(see Fig. 1), and the explicit form of A has yet to be found. All 
that is assumed at this point is that A vanishes at any point of 
Sndy where Nen vanishes, or equivalently where N=. 

The general balance equation postulated reads naturally 
im nonlocal form , 


dean (Ueno, (1) ep (ee enn (2.9) 


(a 
ct 


where V is any region fixed in space. For this to hold for every 
region V not intersecting the interface at time t, it is necessary 
and sufficient that the usual bulk balance equation 


Mec iow) cre (2.10) 
at 


be satisfied at every point not on £,. We note that spatial or 


t 
time derivatives are nowhere meant in the sense of distribution 
theory in this work, so that eq.(2.10) is used only at points not 


lying on x,. It will happen that expressions like Ved, which may 


t 
be undefined on vee appear as parts of integrands in integrals over 
regions containing part of ree Their lack of definition over these 
sets is however clearly immaterial because these sets are of (three- 
dimensional) measure zero. 

Consider a motionless region V intersecting the interface 
De at time t and extending into Me and My as) ti lustrated invrig. Ie 
Our task is now to evaluate the terms of eq.(2.9) for this type of 


Veqione First 1b 1s eastly seenm that 
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where Vv; = VaM, , [al is defined by eq.(2.4), and v- is the normal 


speed of Dy defined by eq.(1.16). To compute d A*(V) we also need 
dt 


the time derivative of the second integral in eq.(2.5). At the end 


of this section we shall prove that it is given by 


ions 
| a> dA = | poea 2 Five dA =) loevea> (2.12) 
t }eedt 
2 (V) 2, (V) L,(3V) 
The meaning of the various terms occuring in this equation are as 


follows (see Fig. 1). z,(V) = Vnrt, is the part of a contained in V 


t 
. The derivative d- is the normal 

Y dt 

time derivative defined by eq.(1.17), H is the mean curvature func- 


and dA is the area element on = 


tion of Dy (defined by eq.(1.13)), and v~ is the normal speed of 


defined by eq.(1.16). L, (av) is the curve bounding the patch z,(V) 


t 


or equivalently is the curve along which 3V meets with rye At each 
point p of this curve, the geometrical factor t is defined by 
y= Net = cotan(.a) (2.13) 
Ney 


where N is the outward unit normal to dV, v is the normalized 


projection of Non Fs , n is the unit normal to r N,n) is 


) a 
the angle between N and n, and dd is the arc-length element on 


Now according to eq.(2.8), the second term in eq.(2.9) 
is given by 
$*(3V) = | jeN dA + | Geta 

vs (2.14) 
“Jt dV + | vedo dV + | (Md -nD+Divj}dA+ | Ad, 


V 
+ 
V. Vi, 2, (3V) L, (dV) 
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Figure 1: The dividing surface a intersecting a three-dimensional 
region V fixed in space, with boundary 9V and outward 
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unit normal vector N. 
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where use has been made of the divergence theorem in three and two 
dimensions. 
Combaningmthise with) eqs. (2.5 jer (ZR 11s C2512)eand 2-6)" 5 


eq.(2.9) may be rewritten as 
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But the first two integrals vanish because of eq.(2.10) and the equa- 


tion reduces to 


figs 
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L, (av) 
When V is so chosen that oV meets Ly at right angle all along the 
curve L, (av), the second integral vanishes and the arbitrariness 


remaining in the choice of V suffices to imply the vanishing of 


the first integrand. This gives the surface balance equation 


- tytas + Divj -o° + Wen-avell = O . (2.16) 
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The function A occurring in eq.(2.8) could a priori depend on any 
property of the surface S defined along the curve Snz, . However 
substituting for V in eq.(2.17) a special region such that part of 


dV coincides with S and the rest meets as at right angle, we get 


| (4 - tavt) ay = 0 
. 
where [ = Snz,. A little thought reveals that the same equation must 
hold for any arbitrarily small portion of T so that A must equal 


ta°vt . Hence eq.(2.8) becomes 


(Ss) = | Jen dA + | Gisv + tie led) (2.18) 


tions. 
A*(V) = | ahh oe Wee alk (2-193 
V Z,(V) 
P*(V) = | o dV + | o° dA (2e20) 
Z,(V) 
SED) = | JN dA + G-0 + tit) aa (oN 
Vv L, (av) 
jae ae et =P Ne (2722) 
Ney 
da + odie ree (at points not on r,) (2523) 
) 
duameachiveas GaDINi meso oen [ciieavell n=an0 (2.24) 
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As far as we know, the contribution tj. - ta’v’ to the line 


flux integrand in eq.(2.21) has not been reported before. 
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Equations (2.23) and (2.24) express the balance of a scalar 
quantity. The general balance equations for a vector quantity may be 
easily deduced from them. Consider a vector-valued extensive attribu- 
te A with bulk and surface densities a, c, a> and o>. Equations (2.19) 
and (2.20) still hold for vector densities. The flux density in the 
bulk will, as in the scalar case (J-N), be a linear but this time vector- 


valued function of N, SO thacethe bulkecontribution to the flux will be 


$(aV) = | M*(N) dA (De 


where M is a field of linear operators and the adjoint has been taken 
for subsequent notational convenience. In eq. (2.21), the term 

tj” - tv'a> will be replaced by Bi eee As in the scalar case (iow) 
the line flux density will be a linear but this time vector-valued 


function of v, Q(v). so that the surface contribution to the flux 


will be 


At each point p « £, Q is a linear map from Ta(24) into E;. Decom- 


p 
posing Qs (u) into normal and tangential parts it follows that at each 


Dae z, there exists a unique vector a € Tay) and a unique 
: S 
linear operator M> : T+(2_,) + T+(2 such that 
in p : “a .) 4 +) 
Qs(a) = (w3)*(a) + Ti n (2.26) 
for every U « Ta(ey) > where the adjoint has been taken for later 


notational convenience. Combining this with eq. (2.25) we see that 


the general expression for the flux in the vector case is 


OMEN. = | M"(N) dA + | C(M)*(3) + Feun +tvtacd da atl) 
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where again the second integral vanishes by definition when Le and 
eV idol noe wmntersect. 1G follows) that ait é@ is a constant vector, the 
flux associated to the scalar attribute @-A will be 


eo (0 n= rac) dA + | [é-(M>)*(p) + @enT*v +1v 6a] dA 
3V L,(3V) 


| M(é)*N dA + | c(M°(@)4eenT)> 0 + tvt@ea] dA 
a L(V) 


where @ = P(é), so that the bulk current density associated to the 
scalar density ea is M(é) and the surface current density asso- 


ciated to é-a a> is M°(é) + @enT . It follows that the bulk balance 


n 


cS Ge =e 
equation for esa 1 


which, since Ve-[M(é)] = (veM)-e + Tr(Mve) = (VveM)-é@ and since é@ is 
arbitrary, implies that 
BEL oa) es OE (2.28) 
ot 


This is the general balance equation for a bulk vector-valued den- 


sity. It follows also that the surface balance equation for the sca- 


d'(€-a°) - 2Hv é@-a> + DiviM°(é)+enT]-e-o>+(n-M(@)-e-av'l = 0 
a (2.29) 
Now we have from eqs.(1.9) and (1.15) 
Div(M°(@)] = (DivM?)-@ + Tr(M> 7é@) 
= (DivM?)-@ + neeTr(MS) 
= @-[DivM> + Tr(SM°) n] 
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But in general 
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UeGrad(e*n) = D-(een) = 


U U 
= = 695(0)e=e0S(e) o0 
where @ = P(e) and use has been made of GS lesen lO mandrel ie)s 
Sor that 

Grad (én) = - S(é) . (2.30) 
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T-Grad (é-n) 
from which 
Diviesn Mra e-f-S(m) > (Divr) nt. 
Substituting these results in eq.(2.29) and invoking the arbitrari- 


a 
ness of e, we get 
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This is the general surface balance equation for vector-valued den- 
sities. To get an intuitive understanding of the terms involving M° 


and T, we note that if ae is a family of simple closed arcs traced 


out on re around a point D; with enclosed area S. such that 5, 0 
aS)e > 0, then 
f = 
lim 1] Q(3) da = Vim 1 | [(M°)"(0) + Ted n] da 
CUS a Vi E20)—S 
€ € Y 


= DivM’ + Tr(SM?)n + (Divr) n — S(T) 
(this last expression being evaluated at Ds 
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F : i or ee cee , 
where in general, given two vectors A,B , A@®B is the linear oper- 
> > aa 
ator defined by A@B(u) = (Bea) A for every U « E3. For convenience 


we bring together all equations pertaining to the vector case: 


(adie | a dV + | ae ay, ea) = | 6 dV +1o6> dA (37) 
V 2 (V) V u4(V) 
— > as 
ey) = | M*(N) dA 4 | ra(s) + rv+254 da (2.33) 

3V L, (av) 

Q = (M°)* + ner (PM? = M>, ne? = 0) (Pou) 
aa+vM-¢ = 6 (2.35) 
ot 


daa Me CHVk = DivM> SeiY(SM?) n+ (Divivin = S(T) = 6° 


(Ou 
cr 


aM (ie vcad = 0. (2.35) 


In concluding this section, we will prove eq.(2.12) which 
is instrumental in many of our derivations. Although related equa- 
tions appear in the literature, as for instance in (Moeckel [1975]), 
we could not find a proof, nor even a statement of the specific 
result embodied in eq.(2.12). 


Let >° be a surface fixed in space and S+ be a piece of 


surface moving about in £°. Given at time dependent function Dp = Foe) 


defined on 2°, it is easy to show that 


r 
d | fdA = | onan J fy a, (2.37) 
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where ve is the outward normal velocity of the boundary curve 


OS, (this velocity is the obvious analogue, for a curve moving on 


a surface, of the normal velocity of a surface moving in space, as 
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defined in Section 1). We shall obtain a result actually slightly 
more general than eq.(2.12) in that we shall find the rate of 


change of the integral 


ne ee [a aA (2.38) 


Und, 


where a iS a function defined on the moving surface r, (that is a 


t 


function a, is defined on each r,s U, iS a possibly changing vo- 


t it 
lume contained within some closed surface ou, the points of which 


move according to some specified velocity distribution Ve and 


U,nx, is the patch of ©, contained in U,. The points in the inte- 


cet t t 


veer One U, need not be attributed any specific velocities as this 


is completely irrelevant to the integral in eq.(2.38). It is also 


intuitively obvious that only the component of V, normal to 3U 


e 
should be involved in computing dl. We will use the normal parame- 
dt 
trization based on re Re » introduced at the end of Section 1. We 


recall that it is defined by 
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where the function R¢(t,p) is such that R+(0,p) = p for every 


De Des and aRt(t,p) is perpendicular to De aca: In order 
to use this St we start by noting the obvious equa- 
ity 
an | a dA | = d_| a(t+t,p) dA (2.39) 
ES eqnmg OTS ot ge coe 


° * 
Let us denote by S. the patch eee es ee at time t+t and by S_ the 


a 


zt tel) 7 soittle iva idioma Bantaeblaortsah eet aie 


sop pntmaend Mefzeogts at yh ay Seana Heri at om ma 
Vat TO Skagway Ate soehior DaAGl I ame Avatiw tantnded at 


‘4% note udeeset S qt gy bet ttange shar ©! ontpioans * 
oadit i ni #tnfoq set . a elniad oD 0° 0 (stey ont ai 
Sint seceatst 20) ae: OF Kidhe wah begufher 4 9c ton bagn yf Yo 10 
Oehd etude Gah. Pha fe hy) be yyaral Bi a? see ei vTedalqnos 2 


se ob Patron y +H Pree an avis oy feet soft ‘wut tydd y Favtrtoant 


a 
<Sratige FeMmiod Ais shu ffiw aw it: gathers if Oey faved 3 binerte 
ee - 
6a. adkese? Yo aha) athtsre Latithe v7 |, BA +s no hbeped wotrosttd 
i) _ oo 


qe barittep sv 3F gale fPeogt 
ey - wl ie a ‘ arly a eo ‘ . «& 
(Giz) Bs 4 a) De cigs gh a 


vhive HOT = (Gh Bh Had wal ei oP “a nedtonst at arsiw 


wh nt 2h, oe Jar etyothhaqnay ebdaen) =! bits wv > q 


Gs, 
~eepe 20al vad otf water es + Ot" a ne! testo iediiea <0 ogu a 


‘aia : 


fit 


2 


r 
a) 


mee i i.) a 


nil = ul ty 3d eres aft 2 


39 


region of De which is mapped into S_ by Rt, that is S* = (2c 


Rromreqs( 22) st follows-tthat 


| a(t+t,p) dA = | a(t+tR°(q)) J, (t4q) Gate 
S Se 
T T 


where uF is the jacobian of the normal parametrization. Combining 


this With eegs<( 2207), (2030). and (2509) yields 


f >| > > 
dI(t) = $—Ta( tet RA(G)) J, (48) dA + | av, dA 
dt St TSO 9S% 
Bi osSee oueba ied ae a (2.40) 
: dt e me ? i 
Und, 3So 


where use has been made of eq.(1.25) together with the fact that 


Se is the identity map on Lee Sovthat os, =U nds. All that remains 


t 
to be done is to compute the normal velocity ue in the second in- 
tegral in eq.(2.40). The curve aS 7 is a closed curve moving, as Tt 
evolves, in the fixed reference configuration ly: lete = Glen) 


be a parametrization of this curve at time t. By the very defini- 


tion of the normal velocity v., we have 


aq = [3 + VV 


OT OT ) 


where V has its usual meaning (outward unit normal to aS tangent 


to ,) and ae = JON mG hale dines, R+(t,q(E,T)) 1S a POUnt of the 
OT 
// 


curve on which ou and ae intersect, and the velocity of this 


ae 
point is 


a R(r,GCE,7))I|_ = BR(r.GCE.t)) + (ROD, (8G(E,7)) 
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where (Ro) is the differential of R (see eq.(1.21)) 
On the other hand, if N is the outward unit normal to ou, at 


the point R(1,.q(E.1)); and V~ is the normal velocity of the surface 


au, ate thats points chien, trom ed ales). 


N+ 9 CR(t,q(E,7))] = vt 
OT 


Using this together with the previous equation evaluated at time 
m= 0sat which Re is simply the identity map on rhe we get 
Vhs vt Neh + vy NS 
: : ae wiped. 
where we have used the geometrically evident equation N- seas GO: 
// 


It follows that 


where now t=Nen is the geometrical factor introduced in eq. 


(2.13) and has nothing to do with the time parameter used above. 


Inserting this expresssion for Ve in eq.(2.40) yields the desired 


result: 
at Al 
d_ | a dA = cs 2Hv-a] dA + | ale een yanigdy G2e41) 
dt paras a2: Now 
ae tget OCU a 
Cert 
where 3(U,nz,) is the curve bounding the patch U,nz,. When U, does 


Chat coe: t 
not depend on time, that is, is a region V fixed in space, Vi o20 


and eq.(2.41), with our previous notations, reduces to eq.(2.12): 


dt 
u,(V) 2, (V) L(V) 


ai 
a dA = | coo 2Hv-a] dA - | tv-a dA. (2542) 
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The three terms on the right-hand side correspond to the three 
causes responsible for the variation of the integral on the left: 
first the time variation of a along a normal trajectory, secondly 
the local stretching or contraction of the surface Z, (see eq.(1.25)) 
and thirdly an edge effect due to addition or removal of strips 

of Z,(V) to the integration domain as its border is delimited by 


varying portions of dV when x, moves. 


t 
The general surface balance equation at the empirical level 
was given perhaps for the first time by Ghez [1966], although his ver- 
sion is not quite correct because he forgot a contribution from the 
bulk density jump. Special versions were given earlier by Scriven 
[1960] and Slattery [1964], while jump conditions in a more restricted 
sense have long been known (Truesdell [1960]). Derivations of the gen- 
eral surface balance equation at the empirical level have also been 
given by Moeckel [1975], Deemer and Slattery [1978], Bedeaux et al 
[1976] and Vodak [1978]. The novel feature of the derivation given in 
this section is that no reference whatsoever is made to tangen- 
tial velocities on the dividing surface, in arriving at eq.(2.24). 
To our knowledge, the result embodied in eq.(2.42) has not been re- 
ported before. From our more general standpoint, we may if we wish 
recover the formulations of other authors, but our formulation can 
also perfectly accomodate the conception of the dividing surface 
as merely a conceptual device, in the spirit of the Gibbs approach 
described in Section 1 of Chapter I. Though of lesser importance, 


the formalism for vector densities embodied in eqs.(2.32-36) also 


appears to be original. 
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3) Mass, momentum and angular momentum balance . 


Again for the purpose of future reference, we shall apply 
the general equations obtained in the previous section to the cases 
of mass, momentum and angular momentum. The system considered will 
be one made of n chemically non-reacting constituants, and devoid of 
internal angular momentum (De Grott and Mazur [1969], p. 305). Here 
too we are careful to introduce the least possible number of assump- 
tions. For example, in considering the balance of linear momentum, 
we could argue as if the dividing surface were material and intro- 
duce a surface stress tensor in terms of surface forces. We could 
then write down at once an equation looking just like Newton's sec- 
ond law. In contrast to this, we merely assume the existence of a 
Q field (see eqs.(2.33) and (2.34)) that gives the surface contri- 
bution to the flow of momentum. This is clearly the least that one 
can assume. The resulting surface momentum balance equation then 
looks rather abstract. To proceed further and obtain an equation 
more reminiscent of Newton's second law, we need some properties 
of the momentum Q field. They could be written down heuristically 
on the basis of an intuitive "material" picture of the dividing 
surface. However it turns out that they follow from the balance of 
angular momentum. We are thus able to show that the intuitively ap- 
pealing form of the surface momentum balance equation follows from 
general physical principles and minimal assumptions which do not 
depend on any specific interpretation of the surface fields in 


terms of the diffuse structure of the interface. 
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The bulk and surface densities of the constituant a of 
the mixture will be denoted by oe and on with corresponding bulk 
and surface current densities J, and ie . The mass balance equations 


are then simply from eqs.(2.23) and (2.24) 


os 
Jka tances 1.0 el) 
a : 
as iPS a 7 > tS 
dey - 2Hv oa Da Jon a irs Ta0e A SE): 


Summation over a = 1,2,...,n gives 


> 
Oper WG a=. 0 (8) 
at 
d'p® - 2Hv'p® + Divg + [Gen - pv'l = 0 (3.4) 
t 
where 
n 2 abe maeiae: n_ 
Cok ROMs pG= pee eGs=meULege=s) Yi, (cS) 
o=1 a=1 @ oa o=1 @ 


We now consider the case of linear momentum. In the bulk 
regions, the momentum density coincides with the total mass current G 
defined in eq.(3.5). The associated balance equation is then 


Gwe F = 0 (3.6) 


where II is the momentum flux density tensor (Landau and Lifshitz 
= 

[1966]) and F the body force per unit volume. The full expression 

for Il in terms of the bulk stress tensor T would be (Samohy1 


[1969], Kehlen and Baranowski [1976]) 


n 
Il = 1+.) 6 W@W 
pepe. er oe 
t > > 
= . Sie 
T+ pwew?e ) Py Wye W, (S52) 


43 


oo 
_ 


v0 rs eer ten wt +4 aettie AgG sone oe rr, 
i pnibinongen > NdYw sabre » Re pe dnieb sf TW wid hm ct 
ana tdbups ‘sonkfad 2esm ont . 1 OS au sagt dtanad Trg Ve ew ‘bn 


ae) 


(45: 0)-bea(ES Rate wort viGmba nedt 


‘- 
FE) wo bw + 08 
\ ak : ne 
; , ' ‘ f ¢ iy . * 
a £} 1 ’ a i tall | : = re v! Ps 
eg i 
; ~ 7 ~- 
| a a) » + 
Gy } : 
#0 
f ei) (Vee i 
1 
‘ re i 
| 
' 7 
~ et T fi 1 ; = ' 4) ‘j elt 
MS TSS eee re ri yd ’ ; s Muswveia 
aus a ee i 1 0349 Sh iA 
’ 
WY 
7 
‘ braeth ite unteeia.) | . Vait ‘ es aA i : 4 wry sy , if or ne t 
CEES 4 " : LA 
oat es bal | Ln gt] aye fo - ae Oo By 604 AGA 
i : | . 
idcmad) a Rivew 7 Yaanet 2agavte dtd ‘sry 
> 
be ( Fee asne ise yaa re ts 
Al 
+ 
i) ’ t an e , t A y _ 
> : \ ae  4) 
—- =h 


= 4:2 u m  & i, 
: 


aa ee : é a4 > See 
where w 1s thes bulk’ Velocity of constituant’ o (J) = PW a wats 
> 
W 


a 
: P > ag > > 
the bulk barycentric velocity (pw = G), andw. = eee Gc In the 


a 


following we shall neglect the last (diffusional) term in eq.(3.7) 
as is done in the second chapter of (De Groot and Mazur [1969]). 
The complete equations could be recovered by considering that the 
last term in eq.(3.7) has been absorbed in - T. So we shall write 
SP <9 ke wie’ woe, (i:)) 
The question is now: what is the surface momentum den- 
sity GP nan obvious candidate would be g + o°v'n, where p> and g 
are given by eq.(3.5). But this requires some justification. We 
will find it in the following general principle. Given a cube C, 
S SoS: 


meee (e 


SS, =O) Wi th 


5 e > > 
that 1S a set of points x = Xo + 7e. 
A Aj 
rj a DOsdie tee 0) denote the total mass current through the 


consisting of all points Xo ec 


Surface S. Nes, S AS se. OC. dil 
Anak AO BERS Sl jy 


(S d 
(0 < X < d) is an intermediate "slice" parallel to these 


and Si; are the two faces of the cube perpendicular to e, 


ONO 
and 37) 
two faces). Then the ge component a of the total amount of momen- 


tum contained in the cube is given by 


Applying this to a cube C the center of which coincides with a 
point p of Z, and no face of which is parallel to Ta(24), it may 
be shown, using eq.(2,18), that the gel component of the total mo- 


mentum contained in C is given by 


P.(C) = | (A) da = é,° | (§+p-vin) dA + | G dV J 


Z,(C) G 
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where Z4(C) is the patch of Z, contained in the cube C. The proof 


is somewhat technical and will not be reproduced here. Hence the 
aS 


surface momentum density G is indeed given by 


pS 


Gedo ven (3.9) 

The surface contribution to the momentum flux follows 
from eqs.(2.33) and (2.34), putting in the Q field appropriate to 
linear momentum, which will be denoted by Q,- The associated M> and 
T defined by eq.(2.34) will be denoted by 1° and Tp so that 

S\x* > = 

(aan lee Thea ne. (3.10) 
From eq.(2.36) the momentum surface balance equation then reads 

dalg ep vn) = chy (Geo vine eDivin ear (Sinan 

dt 

+ (DivT, 

where f is the surface momentum source, which may be interpreted 


)H- S(F,) - f+ cma) - ven = 6 (ii 


as an external surface force per unit area. This equation is hard- 
ly suggestive of a surface version of Newton's second law. In that 
respect the essential element missing is an expression of Tp in 
terms of already defined quantities. The desired expressions will 
follow from the balance of angular momentum. In the following, gi- 


= 


ven A « E3, we shall denote by A* the linear operator defined by 
AX (U) = Ax. The symbol Y will stand for a field of linear oper- 
ators, the one associatea to a given point of space Yo being Yo : 
The bulk density of angular momentum is rx @ and its surface den- 


sity will, according to eq.(3.9), be given by 


iS Bete rx (g + pv n) . (212) 
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The flow of angular momentum through the boundary 3V of a region V 

fixed in space will be given by 

>* (coi sn i f 3s a lias 

6 [= | Gee NA see) Wie Q,(v) vale aeak. (gets) 
OV L, (av) 

It follows that the operator field M of eq.(2.27), which in the 


case of angular momentum will be denoted by Wee may be written as 


Uegece luuidikel) pacet cel Vane 
Sues Ie yr = ies Equation (2.35) then becomes 
mle. (lies wrolenig Vee tae a 


In order to evaluate the middle term, we note that if é is a cons- 

tant vector, then 

lire VeLI(r x @)]= (Vel) -(rx e) CT virx e)] 
S86 Won) Sree) 


Now the antisymmetric part of Il may be written as m for a certain 


vector field 7 and the last term above becomes (because of the anti- 


symmetry of é*) 


Tee ys Iv 6) = =2i-e , 
and it follows that 
ote pS Sees 2a 


Hence the balance equation for angular momentum becomes 
> 


Neale ee ee i Se 


re) 
which in view of the balance of momentum (eq.(3.6)), reduces to 
=~ 


i ete 


Thus the antisymmetric part of Il vanishes and 
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eh lie (3.14) 
In the end the bulk angular momentum balance equation is seen to be 
equivalent to the symmetry of the momentum flux density tensor. 
This is well-known and simply reflects the absence of local torques 
due to the absence of a local intrinsic angular momentum density to 
store up their effects. To derive the consequences of the corres- 
ponding surface balance equation requires a little more work but is 
quite straightforward. The final result is expressed in eqs.(3.24) 
and (3.26) below. First we must find the standard decomposition 
(eq.(2.34)) of a. the surface Q field pertaining to angular momen- 
tum. From eq.(Sel3s)s it follows that 
2S 


Vo. (i) at nel, . (3215) 


Accordingly, for any tangent vector u, 
= -(rxn)-(1j)"@) = (nxt) a 
so that the T surface vector field for angular momentum is 
T= i ee (3.16) 


Likewise from eq.(3.15) it follows that for any tangent vector u, 


j(u) 


POL(G)) =e EP (i) at en) oa 
a Je 
which means that M the M operator for angular momentum, is given 
by 
(3217) 
where use has been made of the general definition of me expressed 
by eq.(2.34) . The above equation is equivalent to 


pee ea toa Tee ee (3.18) 
a Je 
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In order to write down the surface angular momentum ba- 
lance equation, we must first evaluate the terms involving M° and P 
in eq.(2.36) with the particular M and Tr given by eqs.(3.16) and 
(3.18). To evaluate DivM- at a point D of Da we consider its sca- 


lar product with a constant vector @e« Tale If e is the surface 


t): 
vector field defined by @ = P(é), we know from eda ler5 achat 


4): In 


the calculations that immediately follow, the base point Ds at 


Ve = @en S, which at D reduces to (Ve)> = 0 since @« Tale 


which € is tangent to Z,, is omitted but it is understood that 


all the expressions are evaluated at that point. First we have 


+x = 


Diven> Pr’ (8)] - Tri’ Pr ve] 
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But for any tangent vector U, we have from eqs.(1.4), (1.10), 
(Vela rand: (1.15) 


FE [RPP 6 ile Fe Decn Pe aed] 


from which 


and 
Trin? Vop(*#xé)} = 8+ St 
Inserting this in the previous equation for Div (n° Pr] we get 


Divi Pp hole] 9 PL (Divi rd ee ST (Pehle tri Sm enw. =(3.219) 
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Next we need the divergence of the last surface tensor in eq.(3.18). 


On noting that in general 


Div (Uev) = (Divu) V + (Vv)(U) (3.20) 
while 
V(r<n) == (ne + Pr<S) (3221) 


it follows immediately that 
Div(T,@ (rxn)] = (DivT, 


Combining this with eqs.(3.18) and (3.19) gives 


rea nx-p ~ PCrxS(T))] , 
DivM, = P{Fx(Divm-S(F,) + (Divl, + Tr(st?))nJ} 

+ SIP(nx¥) - nxP'p : ie) 
From the general relation Tr(u@V) = U*V and eq.(3.18) we then get 


Tr(SM°) = - Tr(SI°Pr*) + S(T) *(¥xi) (223) 


while using eqs.(3.16) and (3.21) 
Dina Div (1° (nx¥) J 
= nelyxDivil’ 1+ Tr(Si>Pr*) + Trdn*) . 
Combining this with eqs.(3.22) and (3.23) we finally get 
. SS, Sa Sivas Can =p Ay Be = & 
DivM. Pe (oN sone (Divd: ) n - a) = 


rx {Div - Spa (DivT, #Tr(Sm”)1n} z nxT, + Tr(Il 


For briefness let us denote this expression by Q.With the surface 
source term rxf (see eq.(3.11)), the surface balance equation for 
angular momentum,according to eq.(2.36), will read 


dtce(g+ oovin)1 - 2HviLex(g+ psv'n)] + 0 
t 


(oe 


oF + PU(A) - vt aETD = 0. 


However, when it is noted that dat = v'n and use is made of eq. 
dt 


Seo 
tlie ata 
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(3.11), this complicated equation reduces to 


> 


ne Vag a= lp) + Tr(IPn*) n = 0 
which is equivalent to 

Tp = vg (3.24) 
and 

QUI Ol es (3.25) 


Now, any antisymmetric linear operator on Tae may be written 


) 
as an* for some number a so that, denoting by Hie the symmetric 
part of I> and by on* its antisymmetric part, eq.(3.25) becomes 

0 = Tr{(me + an) n°] 

= aTrin*n*] = alr(-I) = - 2a 

SO) that 

(Te me ae (3.26) 
Thus the surface angular momentum balance equation, when combined 
with the surface momentum balance equation, is equivalent to the 
pair of equations (3.24) and (3.26). These equations imply that 
the Q field associated to linear momentum may be rewritten as 
Qe ier vneg 
We now degine a surface tangential velocity field v in terms of 
0° and 9 by writing 

Ge=an ny ae (3-27) 
The above equation then becomes 

Qp Sale ey ov ney ; (3.28) 
Whereas the surface quantities 0° and g may be considered to 


have operational definitions contained in eqs.(2.19) and (2.21) 
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with a> = 0° and j = g, the tangential velocity v has no indepen- 
dent operational definition within the present "minimal" framework. 
It is first introduced in eq.(3.27), as a derived quantity. In an 
intuitive "material" picture of the dividing surface, this equa- 
tion would seem physically obvious and we would tend to consider 
Vas a "primitive observable". We will now see that the surface 
momentum balance equation derived in the minimal theory may be 
cast in just the same form as one would expect on the basis of 
such intuitive interpretations of the dividing surface and of V. 
To this,end a non-tangential velocity distribution V is defined 
on ue by 

ee yen (3.29) 
where v- is the normal velocity of re and V is defined by eq.(3.27). 
When this is inserted in eq.(3.11) and use is made of eq.(3.4) to 


move ne to the left of the time derivative, the following equa- 


tion is obtained 


oC dV + BV +0-S(¥) i + (div + VeGradv') n - vi(Gradv’ + S(¥))I 


Let us denote by D the time derivative along the trajectory Q(t) 
c 


defined by Q(t) e 2, and Q(t) = V¥(Q(t)) for each t, where V is de- 
fined by eqn( 3.29). Ite followse tromeeds..( 1.19 )eGlec0)cande(tec/) 
that the complicated expression multiplying be in the above equa- 


tion is simply equal to DY. Introducing a surface stress tensor 
Dt 
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ee ie oe (3.30) 
the momentum balance equation becomes 


Sp OMMMDAVGe ee TECST ) net ee EET avian Ve a(aiel emm(aesT) 


D 
Dt 
We may if we wish give a very intuitive interpretation of this equa- 
tion. The term on the left may be viewed as mass (per unit area) 
times acceleration. By analogy with bulk continua, the surface 
Stress tensor can be imagined as representing internal surface for- 
ces according to the prescription that the total force is acting on 

a@ patch Of surface S ¢c ry bounded by the closed curve 9S, due to the 


the surrounding "interfacial matter" is given by 


pec lac ede (332) 


This integral is easily transformed into 

F = | [Div eos Tone rues (3.33) 

S 

showing that the first two terms on the right-hand side of eq.(3.31) 
may be interpreted as the resultant per unit area of local surface 
contact forces. The term f may be viewed as external surface force 
per unit area, and the last term as arising form the interaction 
between the interface and the adjoining phases. This informal pic- 
ture iS quite satisfying intuitively. It 1s good to know that it 
may be developed from the minimal assumptions introduced in Sections 
2 and 3, which are themselves independent of such intuitive repre- 
sentations. The crucial step in deducing the intuitively appealing 


eq.(3.31) form the abstract eq.(3.11) was the obtention of eq.(3.24), 
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via the balance of angular momentum. 

At equilibrium, T> reduces to yI, where y is the surface 
or interfacial tension (Aveyard and Vincent [1977]), and from the 
general identity Div(aA) = A(Grada) + aDivA valid for any scalar 
and linear operator fields a and A on x, we get 

DivT? = DivyI = Grady 
On the other hand, it follows from eq.(II.1.13) that Tr(ST°) = 
Tr(Svylje= 24 Ho while ts - pl (U 1s thesadentity on.Es). soi that 
eq.(3.31) reduces to 
Grady + 2yHn + oe (p’-p )nA = 0 
whose tangential and normal projections are 
Grady + f = 0 
and 
p’ = ie pa De 
The second of these relations is the well known Laplace equation, 
which implies that the pressure inside a liquid drop is higher than 
that of the surrounding vapor. 

The balance of linear momentum at the empirical level has 
naturally been considered by all authors who worked in this field. 
We may quote in particular Scriven [1960], Slattery [1964], Ghez 
[1966], Moeckel [1975], Murdoch [1976], Bedeaux, Albano and Mazur 


[1975], Kovac £1976] and Vodaék [1978 a,b]. The derivation given 


in this section distinguishes itself by the minimal amount of assump- 


tions made about the dividing surface. Thus no appeal was made to 
intuitive notions of surface material velocity or of surface con- 


tact forces. As a result, the formalism, which is incidentally com- 
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patible with a "material" conception of the dividing surface, is 
also fully compatible with the Gibbs approach, though neither points 
of view were relied upon in the derivation. This confers a high de- 
gree of generality to the equations. To stress this point we note 
that, if the Gibbs approach were adopted, the surface excess mass 
density 0° would be interpreted as an adsorption (in mass units) and 
as such, with certain choices of the dividing surface , it could 
very well be negative. Yet our momentum balance equations would hold 
true just as well. The price paid for this generality is the slight- 
ly increased complexity of the derivation. 

The surface balance of angular momentum as given for in- 
stance by Scriven [1960] and Slattery [1964] and other authors quoted 
above is used only to establish the symmetry of the surface stress 
tensor. In the new formulation given here, it serves the additional 
purpose of bringing the very general but abstract momentum balance 
equation (eq.(3.11)) into the more familiar form expressed in eq. 
(3.31). We consider it of some interest that this intuitive form 
can be arrived at on the sole basis of a postulated Q field for the 
surface flow of momentum (which is the very least that one can as- 


sume), and the surface balance of angular momentum. 


4) Energy 


Our aim in this chapter has been to give an account of 


the empirical approach to interfacial dynamics that will serve 


as a basis for comparison with the more fundamental theory developed 
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in subsequent chapters. In order to get as clear as possible a pic- 
ture of what can be said at a strictly empirical level, we have 
strived to minimize the number of assumptions and to justify all re- 
sults not by intuitive arguments but by deductions from general 
principles. All our assumptions are contained in eqs.(2.5-2.9). So 
the consequences derived so far must hold in any theory that con- 
tains these assumptions, irrespective of the specific interpretation 
it makes of the surface x or of the surface fields defined on it in 
terms of the diffuse structure of the transition zone. A good example 
of this is the way we arrived at eq.(3.9) for the momentum surface 
excess density, or at eq.(3.31). Using eqs.(3.27) and (3.29), eq. 
(3.9) may be rewritten as 

Gieeeo ave: (4.1) 
Although a priori evident, or so it seems (note that in the Gibbs 
approach, 0° could be negative), we nevertheless showed how this re- 
sult could be derived from a general principle relating momentum flux 
and momentum density (see the argument preceding eq.(3.9)). 

Coming now to the kinetic energy surface excess, nothing 
could seem more natural than setting it equal to ao V2. At the em- 
pirical level, one would indeed be hard pressed to suggest anything 
else. To make the point as clear as possible, the question is this: 
Under the assumption that the total mass M(V) and the total momentum 
P(y) in any region V are given by 

m(v) = [ p dv+ | p® an 
V CV) 


and 
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P(V) - | ow av + | ofan, 
V u4(V) 

can we argue from some general principle that the total macroscopic 
kinetic energy (neglecting diffusion kinetic energy if present) in 
V, K(V), is given by 

KGU ) a= | Tow> dV + Boy? el es (232) 

V zr V) 

Surprisingly the answer is no. This simply cannot be proved and must 
be considered as an additional assumption. As we shall see in Chap- 
ter V, when mass and momentum densities are considered more funda- 
mentally as continuously varying field in the transition zone, and 
surface excess densities are given explicit definitions in a gen- 
eralized Gibbs approach, the above expression for the kinetic energy 
surface excess is found to hold only when certain residual terms 
can be neglected, even in a context where eqs.(2.19-24) are fully 
justified. In contrast to this, the general theory will confirm 
that eq.(4.2) holds true in such a context. This possible complica- 
tion with the kinetic energy is just a particular instance of a gen- 
eral situation: when we construct some bulk density a by a nonlinear 
expression a = f(a, ..+-.a,) involving other bulk densities Apoeees 
a> jt cannot be said that the surface excess density a> associated 
to a is necessarily equal to the same expression wherein the bulk 
densities A, o+++5a, are replaced by their associated surface ex- 
cess densities anoepeea ye (In the case of kinetic energy, the bulk 


> 
densities involved are mass and momentum densities p and G and 
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the nonlinear expression is p7! GG). This complication, to which 
we shall return in Chapter V, does not seem to have been explicit- 
ly recognized in the literature. 


The balance equation for ee follows naturally from 
a 


the momentum balance equation. From the definition of the 3 time 


derivative given after eq.(3.29) and from eq.(1.19), it follows 
that 


Dpo® = d*p? + V+ Gradp® 
Dt dt 


which, combined with eq.(3.4), gives 


Ss Dy 


D 0° V2 - 2Hytp°Vv2 + V20°Divv + V2EG-n - ov I x 
pe (2) pe 2) Be a 


2 t 2 Z 2 
Taking the dot product of eq.(3.31) with V and replacing the left- 
hand side by the above expression for p°D (V2) yields the desired 


Dt'Z 
result: 


d’(p°V2) - 2Hv'p°V2 + Dive pev2¥-T°(¥)] + Trl T°(WW- v's) 1 
2 2 


a 2 


i) Aneta Sonia eum 0 - (4,3) 


Comparing this result with eq.(2.24) enables us to identify the 


surface current density and source term associated with SP Ge 
2 


Dp Veve= To (y)mandees tr ET wnvenes) le f-V. There is no question 


that eq.(4.3) is the correct balance equation for the quantity 


o V2 . But if this quantity is indeed the kinetic energy surface 


excess, its balance equation must also follow from the general 


ey) 


a: 


: ai ~ . . ~: 
ae 7 — 7 a heh 


gina fe ‘atthe it on 
swat tor tt REE Tee 


- pw 


i 


“Sob sah E).pa ribtw tentdimes.; 


OR «LUN apes ~ Cb 4 tba 


- pst ait priisa laa t) wz ‘ ay 4Tt pa Ww fobeng 25 on? ve 
dathesi wry! cbt aty aye “Shot seelaes avere, Bay yo aie 


(gevi<d7} Tt dar 5 Ur Del » ~ - a 


(6.5) y 3 eye ae $oselich tsi ier ’ % torr 


“f 


sia i la 03. Be qatcens (ORS), papiiiy Mart ania cd 
25 » Sa bane soit ae me borne! hog Haas nie 


7 @ 


AdWzsuc 67 4) watt vy ' are aa es ieee a << 


et teeHs: vals eat poannt-s baler 3 3: Pia Bree te. BS | 


he epithe quiets ctzent hl ait, buesef af pod ae Se 
hain ot iota fe im it 


58 


surface balance equation (2.24), feeding in the appropriate bulk 
kinetic energy terms for the jump [Jen - av" I . The bulk current 


> 
density J associated to kinetic energy is ow2w - T(w) while av~ 
2 
becomes ow2v". From the above identification of surface current 
2 
density and surface source term for ae. there follows an alter- 
2 


native surface balance equation for ove : 


tel - 2Hvp V2 + Div Co V2¥-T°(¥)] + TrLT*(W- v'S)] 
t 


It is compatible with eq.(4.3) only if 


[-V-T(n) + (Gen- pv')(Vew-V2)2= Clow2w-T(w)Jen - pw2vtd , 
Ze 2 2 


which may be cast into the simpler form 


2 xe = 
neQ(T - ow U)(w, ) DI = 0 (4.5) 
where 
= nh 
yew - (4.6) 


and U is the identity operator on E3. This conclusion seems 


inescapable inasmuch as we accept first oie ee the correct ex- 
2 


pression for the kinetic energy surface excess and secondly the 


surface current density and source terms oVv2¥ - Ty) and 
2 


Set (iy ES aes f-V as those appropriate for the kinetic 
energy. By virtue of the generality of our formulation, it follows 


that this state of affairs obtains no matter what specific inter- 
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pretation is given to the dividing surface and to the surface 
fields. But it has been overlooked in all the studies that we 
know of where the surface balance of kinetic energy was consid- 


ered. 


ata AE 


> 3. OS 
*“n=Wen = V on 2} 


Equation (4.5) is satisfied if w ‘ 


and either n or both W, and Wi, are eigenvectors of T randevs 

on rye This is a sufficient but of course not necessary condi- 
tion. Another sufficient but not necessary condition is that 
Ww o=w = Von re We shall get better insight into this matter 
when we consider it again in Chapter V from the point of view 
of the general theory developed in the next two chapters. 

The balance of kinetic energy is of course an inter- 
mediate step towards the balance of internal energy which in 
turn leads to the balance of entropy and to the expression 
for the entropy production, which is the starting point for the 
derivation of linear constitutive relations in phenomenological 
non-equilibrium thermodynamics. Since the results thereby ob- 
tained are rather involved, and since we shall not need them 
in subsequent developments, the reader is refered to the litera- 
ture quoted in Section 2 of Chapter I, particularly to the works 
of Bedeaux et al [1975], Kovac [1978] and Vodak [1978 b]. 

In bulk phases, the balance equation for combined 
kinetic and internal energy e = u + dpw*, ignoring external 
fields for simplicity... 1s 


ee | ae Pe le (4.7) 
aE q 
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On substracting from this the balance equation for kinetic energy 


3 (ow?) + V*(ow2w-T(w)) + Tr(TVw) = 0, 
at 2 2 


there follows the balance equation for internal energy 


ie Ve (uw+ JQ) = Ty ( Ty) =) Ole (4.8) 


Q 


In writing down the corresponding surface balance equations, all 
authors proceed by analogy. In our notations, this means that the 
balance equations for the energy surface excess is written as 


dte® - 2Hvte® Dey le ance ] (4.9) 
dt 


+ Tew- TOW) + d0) i Prey 1-20 
which, combined with eq. (4.3), yields the following surface ba- 
lance equation for internal energy: 


dr a ene = is aS eG esi (4.10) 


dt q 


Now it should be stressed that eq.(4.9) is obtained on the basis 
of analogy and analogy alone. Once we have identified e° as the 
surface energy excess and jg as the surface excess heat current 
density, there is in principle no freedom left in the choice of 
the surface excess energy current density. That it may be written 
as ev-T°(V)+Jj. is an additional assumption. We emphasize that 
contrary to eq.(4.1), we have not deduced this expression from a 
general principle. If the surface 2 were a model of a free liquid 
film, we could in fact justify this choice on the basis of a work 
principle. But in the present case we, at least, cannot think 


of any way of doing this. 
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CHAPTER IIT 


THE DEVIATION APPROACH 


1) Towards a more general theory 


In the previous chapter, we have aimed at presenting the 
phenomenological approach by starting with a set of minimal assump- 
tions and trying to develop the rest either as direct consequences 
of these assumptions or of these assumptions coupled with well- 
established general principles. This proved possible except for 
the eneray-related equations which required additional special 
assumptions whose sole basis was their plausibility. The approach 
followed throughout Chapter II was fully phenomenological in the 
sense that no reference whatsoever was made about the interpre- 
tation of the surface fields in terms of the microscopic diffuse 
structure of the interface. 

In Chapter I, we saw that the Gibbs approach provided 
a beautifully simple and yet exact rational basis for explicitely 
interpreting the equilibrium surface excess densities. In view of 
its attractiveness, it is worth studying the possible extension 
of this method into the realm of non-equilibrium. It is to this 
matter that the rest of this work is devoted, while the first 
part will provide us with a background against which the forth- 
coming results may be conveniently compared. We shall now describe 


in broad terms what these results are all about. 
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The object of our considerations is a non-equilibrium 
n-component two-phase fluid system. We now explicitely recognize 
the three-dimensional character of the transition zone, and let 
the system be described by densities, current densities and source 
densities, typical ones being denoted by ae J* and o, that vary 
continously across the interface. With the help of these densities 
one may compute the total amount A-(V) of attribute a” in a region 
fixed in space V, its flow through 3V, o (3V), and its rate of pro- 


duction in V, P (V). These quantities obey the balance equation 


eA C UNE onto) sepa (Viera se (cial) 


oO 
ct 


As in equilibrium, the first step consists in introducing a divid- 
ing surface located somewhere in the transition zone, and in extra- 
polating the bulk densities in both phases up to this dividing sur- 
face. This means that they are subject to the usual balance equa- 
tions and that they obey the bulk constitutive relations right up 
to the dividing surface. We call these densities, which in general 
differ from the exact densities in the transition zone and may pre- 
sent jump discontinuities at the dividing surface, the schematic 
densities. The fields equal to the difference between exact and 
bulk densities, we call deviations. They vanish outside the transi- 
tion zone. 

We next come to the second step which involves two parts. 
The first is to introduce surface excess quantities that will com- 
pensate for the errors following the use of schematic densities in 


computing the total amount of an attribute in a region V encompass- 
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ing part of the dividing surface. The exact relations obtained in 
this context differ from those used in the empirical approach but 
complete agreement obtains in a certain well defined zeroth order 
approximation. The second part concerns the dynamics. Equation (1.1) 
is satisfied for any region V not intersecting the transition zone 
provided that the bulk densities obey the usual differential balance 
equations. In the transition zone, the exact densities are the sum 
of the schematic densities and of the deviations. We will obtain 
a necessary and sufficient condition in order that eq.(1.1) shall 
also hold for any region V containing part of the dividing surface. 
This exact relation involving the bulk field jumps and appropriate 
Surface moments of the deviations is goxmally identical with the 
surface balance equation used in the empirical approach. We shall 
explain why this formal identity is remarkable and see how the dif- 
ferences that remain in the detailed interpretation of the surface 
fields disappear in a suitable zeroth order approximation. 

A basic premise of our theory is that the detailed struc- 
ture of the transition zone may be described by exact densities au 
J* and o” obeying the relation 


Sig 20 (12) 


everywhere including inside the interface. This type of relation 
is one basic ingredient in the phenomenological treatment of bulk 
fluids. As is well-known (Kreuzer and Beamish [1977]), such rela- 


tions may be derived from microscopic theory using reduced distri- 
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bution functions and the first equation of the BBGKY hierarchy. In 
the present context, the status of eq.(1.2) is somewhat different. 
For the phenomenological approach to interfacial dynamics makes no 
use of exact densities and consequently of eq.(1.2) in the transi- 
tion zone. It rather postulates surface excess densities and derives 
a sutface balance equation (eq.(I1.2.24)). So our justification for 
using eq.(1.2) in the interface rests directly on the microscopic 
theory. Specifically, what we assume is that the microscopic theory 
of reduced distribution functions can provide (probalistic) ex- 
pressions for mass, momentum and energy densities and current den- 
Sities, valid in the transition zone, and satisfying balance equa- 
tions like eq.(1.2). The use of reduced distribution functions in 
the interface is in fact common practice in the statistical mecha- 
nics of the equilibrium interface (Navascues [1979]), and there is 
no need here to go through an explicit microscopic derivation of 
the balance equations, because the calculations are essentially 

the same as for bulk phases. The only difference is that lineariza- 
tion of the pair distribution dependence on the center of mass 

of two molecules (Green [1969],p.36) cannot be made in the transi- 
tion zone. But that does not prevent one from arriving at the ba- 
lance equations (see for instance Gray [1968],p.514-520). In a 
paper mainly devoted to equilibrium questions, Carey, Scriven 

and Davis [1978] have given such a derivation for the case of mo- 
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2) Density integrals, currents and surface normal moments 


In this section, we give the explicit construction of 
the surface quantities in terms of which the theory is formulated. 
We introduce "surface normal moments" as convenient building 
blocks. These are used to define surface excess densities which 
are the exact analogues of the "adsorptions" defined in the equi- 
librium Gibbs approach (Defay et al. [19661]), and surface current 
densities which are not present in the equilibrium theory. The 
rigourous expressions for the total amount of an attribute ina 
region V or its flow through the boundary 93V are not identical 
with those found in the empirical approach (eqs.(2.19)and (2.21)), 
but complete agreement is obtained in a suitable zeroth order 
approximation that will be discussed later. However, as we shal] 
see in section 3, a rigorous surface balance equation is arrived 
at which is formally undistinguishable from that obtained in the 
empirical approach. This is quite remarkable. 

Many notations are carried over from Section 2 of Chap- 
ter II. The dividing surface at time t is denoted by Dye 1S) Una 
normal by n and the regions My and My have the same definitions 
as in the previous context. The densities a, o and ae defined 
everywhere in the material in terms of the extrapolated bulk den- 
Seas le UU) ea CGO he Te My and a(y) = a (¥) if fe Heo 
(with similar definitions for o and wh. will be called schematic 


densities. Their jumps at the dividing surface are again defined 


by eq.(II.2.4). The name deviation will be given to the functions 
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defined everywhere in the material by 


aLe=na “ciNeeOn ec cand A) a eee C25) 
It follows from the definitions that 
feo ae SG ei (2.2) 
ot 
everywhere in the material, while 
egy ce=e Oh and acao ar eee ro (23) 


ot ot 
in the interior of M, and of My . The same comments that were made 
after eq.(II.2.10) about the meaning of derivatives still hold. 

A convenient and physically reasonable hypothesis about 
the deviations a',o' and d' is that they vanish everywhere except 
in a very thin layer enveloping the dividing surface. Strictly 
Speaking this hypothesis might not be verified in an exact solu- 


LION sa Orminstancess | fOr some D eo ee thes une ti Onee a'(p+tn(p)) 


t? 
happened to be analytic on the interval (0,°), it could not vanish 
on some subinterval (Zo,°) without vanishing on (0,91 as well. 
Granted that such a possibility definitely exists, it is never- 
theless certainly still safe to assume on physical grounds that, 
for c outside some very short interval (-6,6), a'(p+ttn(p)) is 
completely negligible. Should that not be the case, this whole 
scheme would be of no interest in the first place: For the sake 

of technical simplicity then, the convenient hypothesis of strict 
vanishing of the deviations at a very short distance from the divi- 
ding surface will be retained. 
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Given a surface © with unit normal vector n, and a num- 
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ber c, the surface parallel to XY at distance zc, ae » 1S by defini- 
tion the set of points p+cn(p), where p ranges over = (O'Neill 
[1970]). The above hypothesis on the deviations may be reformulated 
thus: it is assumed that there exists a small distance 6 such that 
outside the transition region I, (I for interface) sandwiched bet- 
ween the surfaces ee and oe parallel to vee the three deviations 
asic and J' vanish. The same value of 6 may be assumed to work 
throughout any time interval of interest. 

Naturally the dividing surface is chosen so that it bends 
sensibly only on a macroscopic scale, that is the normal curvature 
Kp () in any direction 0 (see eq.(II.1.12)) is negligible compared 
to 1/6. Its location is further restricted by the requirement that 
the variations of the exact densities along the dividing surface 
must be small compared to those along the normal: variations along 
the dividing surface are on a spatially slow, that is macroscopic, 
scale. This is tantamount to the requirement that the direction of 
n(p) be close to those of the gradients of the exact densities at D. 

The following symbols for density integrals will be adop- 
ted: 


A(V) = | aur, Ae = | a dV ANC a= | at ad) See2 2) 
V V V 
with similar obvious definitions for the integrals of the source 
terms P(V), P (V) and P'(V). Likewise we shall write 


o(aV) = [3 dA , o*(aV) = [ o-i dA , o'(aV) = [oi cA enc25) 
oV OV OV 
To avoid repeated circumlocutions, we shall further agree 
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on a set of notations used throughout the rest of this work. Some 
of the objects defined are already illustrated in Fig.1 of Section 
2 of Chapter II, while the others are illustrated in Fig.2 . Consi- 


der a 4ixed region V intersecting ©, and extending beyond the tran- 


ie 
sition region I, (defined above), into My and M, . The outward unit 


nermal to aV is denoted by Ns the patch of Le contained in V, by 


r,(V); the boundary Cz, (V)] of this patch, which may be described 


| 


as the curve along which dV meets with £,, by L, (av); and the strip 


ie 
of the surface dV lying between the parallel surfaces ry° and Ses by 
A, (av). The "belt" consisting of all points p+etn(p) with p « L, (av) 


and z « [-6,6], is denoted by g, (av), and the part of I, circum- 
scribed by and includind this belt,by I,(v). The unit normal to 

dt c a é Feet : 

B,(3V), pointing away from I,(V), is denoted by N*. The region boun- 
6 =) 
t and Le A 
6, (av), is denoted by X, (av); its cross-section looks like the union 


ded on top and bottom by = and on the sides by A, (dV) and 


of the top and bottom surfaces comprised between the bars of the 


Jeuceue nn cTUNeCtCION Clas) defined. skirst In thewinterlOr.0te s(oU.). 


15 
by e(r) = -1 if F is inside I,(V), and e(r) = 1 if F is outside 
t 


I,(V); it is then extended to the boundary of X,(9V), deprived of 


ic 


L,(0V), by continuity. (The function © usually has a jump when r 


+ 
goes from the + to the - part of X, (dV) .) The outward unit normal 


to the curve L, (av) is denoted by v; this vector is tangent to Lye 


The direction of T, unit tangent vector to L,(3V) is chosen so that 


t 


dxT = n. Since the curve L, (av) laese ano. .1t follows that N is 


perpendicular to T, so that N lies in the plane spanned by n and \. 
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Figure 2: A portion of the region X,(3V). 
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Finally the geometric factor t is defined on L,(dV) as in equation 


i 
(1is2213)%) Fromethesdefinitions,.it issobvious that 


Aa (CU)eeoAt eating dV eed. das (2.6) 
1, (V) X, CoV) 


The integral over I,(V) may be done by iteration: first along the 


6 


normal and then over the surface z,(V). Let r(&1,&2) be a parametri- 


zation of 5, or part thereof, such that ne(dr x dr ) > O . The neigh- 
: aE? 


boring three-dimensional space may be parametrized by 
GE Ae ye) tied erate ne 
0 


> 
In terms of U; = x= and U2 = ee the area element on ae may be ex- 


pressed as dA = ne(U,xU2) d&'d=?. On the other hand, it follows from 


eq.(I1.1.10) that 3R =u. -zS(U_),while 9R=n, so that the jacobian 
J€ oC 
Je ak, (aR x ak is given by 
aw ae an 


J = We{iyxU2-cLU1xS(U2) -U2 S(O, ) 1452S (Hy )xS (U2 )} 
= = (Uy xo) (4 2cH + <°K), 
where use has been made of eq.(II.1.13), and the volume element is 
GV Gawd de de der =\4-—2cH+ ek) dada Cane 
where dA is the area element on rye This is the requisite factori- 
zation. 

At this stage it is convenient to introduce quantities 
that will be called normal moments. Given a function f defined 
throughout the material at each instant t, the normal moment of or- 
der =m of f, “te » 1S the function defined on Dy by 


6 
<f> (B) = | o™eCBacicBy) de (2.8) 
=f 
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where m = 0,1,2,... 1S any positive integer and D e Z,. The normal 


t 
moments of the deviations a', o', J (defined by eq.(2.1)) turn out 
to be natural building blocks for the theory. They will be attri- 
buted special symbols 
> aa 
Ohba de te Otel Sie ww ed Uae (2.9) 
The vector J, will be decomposed into its normal and tangential 


Dac Sars 


“te = al ee Sue aed 
n 


Veneer We ins =O eslagre W240) 


We emphasize that Cry Che and J, are functions defined on ty 


only. Many equations will be simplified by introducing the fol- 


lowing functions also defined on ae 
Soy = 2Ha, + Kap (Zettel 


Ole = Og = 2Ho ate Ko»5 C2) 
For a reason that will be clear shortly, they will be called 
Aurface densities. 


Using eqs.(2.7) and (2.9), it follows immediately that 


r 
avd \ a (a -2Ha,+Kaz) dA = | a. dA 


J 
1 (V) LCV) L(V) 
so that eq.(2.6) may be rewritten as 
PF 
OUD © NGG =: | adi + | oa’ dv . (2.13) 
r40V) X, CoV) 
An identical equation obviously holds for the internal production: 
DAG = Tay | o. dA+ | Go' dv. (2.14) 
zr CV) X, CoV) 


Thus the total amount of an attribute in V is seen to be 
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the sum of a bulk contribution from the schematic density, of a sur- 


face excess and of a fringe effect. The latter, which is represented 


by the last integral in eq.(2.13), depends on the shape of A, (0V); 
it vanishes when A, (av) = 8, (3V) 


An equation similar to eq.(2.13) will now be obtained for 
the flux 0” (9V). The starting point is the obvious relation 

@°(3V) = (av) + (av) . (2.15) 
Since the deviation J! vanishes everywhere on dV except on the 


strip A, (av), the last term may be rewritten as 


[as as 
o'(0V) = J'-N dA 
A, CoV) 
In order to extract a line flux through L, (ov), this equation is 


rewritten as 


DUCA UO t SN adh ell Neda eet gd °-NTocAy |. (2.16) 
dh ala 
B,(3V) A, (av) (av) 


The first term on the right can be expressed in terms of the first 
two normal moments of J'. Let r(x) represent the arc-length para- 


(av), with orientation such that dr = T, where T 
da 


has been defined previously. An obvious parametrization of 8, (av) 


metrization of Ly 


is then 
(Asc) = (A) + cna) . (aed) 
It follows from eq.(II.1.10) that 9p = T - cS(T), while ap = n 
oA oS 
so that 
a) ee 1 BS os 
N* = | Oe ¢ z | oy c 
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Es) 


= | LT-cS(T) 1xn| ; PG ie hi 
or 
(eae, Usa aT (2.18) 
|T-2S(T) | 
Using dA = |dpxap| dadz = |T-cS(T)| dadz together with eqs.(2.9) 
dA dF 


and (2.10), we then get 
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| J'-N* dA = | da | de |T-cS(T)| J'-t-n » LES 
(av) L.(avy -6 een) 
= | io dee [5 jietnxS(T)] da . (2.19) 
L,( dV) L,( dV) 
Now if U is any vector tangent to zr, at a point of L, (av), it fol- 
lows, using the self-adjointness of S, that 
U+EnxS(T)] = S(T)-(Gxn) = T+S(axn) 
= (xT) +CAxS(Uxn) J = -0-0fxS(OxA) J 


It is easily seen that the linear operator Gms -nxS(Uxn) has 
the same eigenvectors as S but with eigenvalues permuted and multi- 
plied by -1. Hence, from eq.(II.1.13), 

So sees ei (2.20) 
where I is the identity operator. The above equation may then be 
rewritten as 

WO = eS) Se (aes 
The operator S will occur frequently in the following developments. 
Its physical dimension is that of a curvature. We note that, since 
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so that 
S Sime) Bee Ole 02822) 
Combining eqs.(2.19),(2.21) gives 
oe dh | eee eh (20123) 
g,( av) L(V) 
Let us now return to eq.(2.16). By considering separately the parts 


of X, (av) contained in My and M, , dividing each if necessary into 


ite ) 
a finite number of regions throughout which the function © is cons- 
tant, and remembering that the deviations vanish on ie and Bs ie 


is easily shown that 


= So ss f > 

NY rahe Se | Aah QO Ved' dV 
ale 

A,(V) B,( av) X, (dV) 


Inserting this together with eqs.(2.16) and (2.23), into eq.(2.15) 
yields 
Do (3V) = OV) ee Igoe San) leo des loved' a. (2a2 a) 
L,CaVv) X, CoV) 
This equation is an intermediate result to be used in deriving the 
surface balance equation in the next section. In preparation for 
this, one last calculation will be made. From a well-known theorem 


about integrals over moving volumes (Prager [1973],p.75), it fol- 


lows that 
d {oatw = | 9 22 av + | oa'w dA (2.25) 
t 
X, (a) K, (av) aX, (3V)] 


where w is the normal velocity of the surface a[X, (av) 1, counted as 


positive when the enclosed volume expands. The surface a[X, (av) ] 
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is made up of A,(9V) and 6, (3V) on its sides, and of parts of Be and 


—6 
t 


A, (av) w vanishes because V/ is fixed in space. This implies that 


1B 


2, On its top and bottom. On the latter parts a' vanishes, while on 


Oa'w dA = Oa'w dA . (2226) 
aLX,(V) J 8, av) 


In order to express this integral in terms of normal moments, we 
must first find an explicit expression for the normal velocity w of 


the "belt" 6, (av). Intuitively one would expect to see an 
dt 


( = -Gradv~ ) coming in. But the final result is not at all obvious. 
Let ty be a given instant at which the integral of eq.(2.26) is to 


be evaluated and introduce a parametrization x(A,n) of the strip 


Ay (3V) such that r(A) = x(A,0) is the arc-length parametrization 
> 
of L, (aV) and T+5% (2,0) = 0 for every X. When X is held fixed 
0 é 


and n varies, x(X,n) traces out on 3V a curve meeting L, (aV) at 


to 
right angle. This curve will be called the i\-curve. Now as t progres- 
ses in a small interval (to-e.to+e), the curve L, (ov) moves on av. 
For each value of X, the point of intersection of L, (av) with the 
A-curve will be denoted by O69 At any fixed t in (t)-e,to+te), 
the map > Q, (t) is a parametrization of L,(V), from which we may 
construct a parametrization of the "belt" 8, (av) 

canoe NCO Rem CNGS 9 Be (2.27) 
To proceed further, we need the velocity at t = ty of the point 
Q(t) (with fixed X}). To this end, we momentarily abstract ourselves 


from our present calculation and consider the following question:gi- 


ven a moving surface oe and a motionkess curve T always piercing the 
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moving surface ©, at only one point, what is the velocity of this 


UG 
point of intersection ? To answer this question, let us introduce 
the arc-length parametrization r(s) of fr and assume that the moving 


Surface 5, at time t is represented locally as the set of points y 


t 
Such that 
AGU ie = eae 


of 
It follows immediately from this that v= pales 
|VF | 


Let us denote by s(t) the [-arc-length parameter associated to 
the point of intersection of I and ae Since r(s(t)) is on De 


at each time t, it follows that f(7(s(t)) ,t) = 0 for each t. from 


which 
Sdpour. oF = 0 
ds ot 
i vo 
which, together with the above expression for v gives $ = = = 
nedr 


ds 


e 
fate. 
12) 


Coming back to our problem, the curve IT is now the }- 
> 
curve and the moving point of intersection is Q(t). We want its 


velocity at time t,. The vector dr of the general situation above 


ds 
is now clearly Nx T and it follows from the above equation that 
dO, (ty ) : (oe Nie 
oe Hie (NT) 


> > a Ree 


qe 
Now we have N = (Nef) n + (N°0)O from which NxT = -(Nen)o + (Ned) 
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and finally 


dQ, (to) - v'(n-ty) (2.28) 
dt 


where « ds defined by eq.( 1122.13). Fromveqase(11.1.527).(2.27) 
and (2.28) it follows that 


dO to sdet u=V (n-tv) 2 cl=Gradvo¢ rveSiv) i. (2.29) 


Q 


It is evident that, when t = t,, the parametrization 6 coincides 
with that defined by eq.(2.17), so that the outward unit normal on 
the part Bp, (3Y) OF aX, (aV)] is -ON*, where N* is given by 
eq.(2.18). It follows that w may be written as 


w = -ON*+ 90(ty 3A5Z) 
at 


- Q Crv+(1-27H+02K) + 0e(c1+22S)(Gradv-)] , 
|T - cS(T) | 


where use has been made of eqs.(2.29),(2.21) and (II.1.13). Inser- 
ting this in eq.(2.26), while remembering that dA = |T-S(T)|dAdz on 
8, (2V) » and combining the resulting equation with eq.(2.25) 

yields the desired result: 


dt 
X,(2V) X, (OV) L,(3V) (2230) 


3) The first surface balance equation 


As the title implies, the balance equation derived in 


this section is not the only one of its kind. However it is the 


d | @a' dV = [032 dV + | {tv'a, + UeLaiGradv” + aoS(Gradv~) ]}dd. 
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most important one because it connects the jumps of the schematic 
densities with the surface excess densities. The balance equations 
derived in Chapter IV do not involve the schematic densities and 
they apply to. higher order normal moments. It is remarkable that 
the first surface balance equation, eventhough it is obtained 
without making any approximation, is formally identical with eq. 
(II1.2.16) obtained in the “empirical” approach. The differences 
that remain concern the physical interpretation of the surface 
densities entering into the equations; they disappear in a suit- 


able zeroth order approximation which will later be discussed at 


length. 

The first of eqs.(2.3) guarantees that the balance equa- 
tion 

a) + o"(av) - P*(v) = 0 (Sen 


holds for any region V fixed in space, not intersecting the tran- 


sition region I,. The first surface balance equation is essen- 


t 
tially a boundary condition equivalent to the requirement that 
eq.(3.1) shall also hold for any region V fixed in space and en- 
compassing part of the dividing surface. 


The first term in eq.(3.1) may be computed with the help 


of eqs.(2.13) and (2.30), together with the two relations 


+ Sy 
d | adv = | el zip | “8 aves | Tal v dA 
dt - i r (Vv 
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where _ Stands for the normal time derivative defined by 


ean tial and 7a1Sedef ined bysecaa(1 lo 2.4o) = ihe tinstuot tiese 
two equations is obvious and the second one follows from the gen- 
eral identity expressed by eq.(II.2.42). Using eqs.(2.3) and 


(2.23), the flux o*(3V) may be written as 


A | enol oe | aoa. | {(lden + Div£3o+S(31) 1} dA 


+ = 
Ue ve r(V) 
Ja’ 
re =16, )e dV 

X, CoV) 


Substituting the last three equations together with eqs.(2.14) 


ande(2.30)i intoveq: (341) yields 


{=—s - 2Hv"a. + Div Hats (ga arcrad Ves (ascrvad vey d 
Da) 
- o, + MJ-n-av" I} dA = 0 


where we have set IJ Ien Sec = [den-av- . In view of eq.(2.3) 


and of the arbitrariness of V, this equation is equivalent to 


déa. - 2Hv"a. DIN (ince jaan Grad ¥-4.5 (ar Gradve ll 
dit 
-o. + ieneave Tero (3.3) 


which is the first surface balance equation. Upon writing 
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it takes the more compact form 


drad--echvealee Diva soul ie nea vem =n0l. (205) 
dt” S S 


We note that this equation is formally identical with the sur- 
face balance equation obtained in the empirical approach discussed 
in Chapter II (eq.(1I.2.24)). That a balance equation of exactly 
the same form as that found in the empirical approach should hold 
rigourously in the present theory is quite remarkable. For it must 
be remembered that eqs.(1II.2.19-22) which define the surface den- 
Sities in the empirical approach do not rigourously hold in our 
theory; they are replaced by the more complicated eqs.(2.13), 
(2.14) and (2.24). The crucial step in the derivation was the de- 
composition effected in eq.(2.16), leading to eq.(2.24), the last 
term of which had just the right form to cancel the fringe inte- 
grals (integrals over X,(dV)) appearing in eqs.(2.14) and (2.30). 
This welcomed formal identity of the surface balance 
equations does not mean however that everything is the same in 
both descriptions. For one thing, if we look at the surface cur- 
rent density j appearing in the surface balance equation of Chap- 
ter II (eq.(I1I.2.24)), we note that it is the same one that is 
used in computing the surface contribution to the flux ®*(aV) in 
eq.(II.2.21). But in the present theory, things are not so simple. 
The surface current density j occuring in eq.(3.5) contains four 
terms explicited in eq.(3.4), and these do not at all appear in 


eq.(2.24) which moreover has a more complicated structure than 
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eq.(I1.2.21). To make a definite comparison, we consider a region 
V, fixed in space, such that at some instant tA, (av) = 6, (av) 
(in general, V being motionless, this relation cannot be preserved 
in the course of time). From eq.(2.24) we have, for such a parti- 
cular negion at that particular time t 

So (eUe= O(oV) + | ijn poly a) 

L, CoV) 

without any approximation. On the other hand, the surface current 
density j which occurs in the first surface balance equation 
(eq2(3.5)) is, according to eq.(3.4) 5 given by 

FS Jo: + S (jx) * aiGradv- + StasGradv-) , 
which differs from the above integrand by the term eaGiad yaa: 
S(a,Gradv-). This is in marked contrast with the empirical approach 
wherein the surface current density used in computing the flux 
and that occurring in the balance equation are one and the same. 
The physical interpretation of the term a,Grad v~ ~ S(acGradva) 
is easily obtained. At a given instant t,, let us trace a simple 
closed curve i on ra Next let this curve be carried along 
with the dividing surface in such a way that the velocity of each 
of its points be always perpendicular to the dividing surface. 
A new closed curve The is thus obtained on each rye 
the "belt" constructed along the curve Ir 


Let ewe) be 


t in the same way that 


6+ (3) was in section 2, and let I(r) be the portion of the tran- 


sition zone 1 circumscribed by this "belt". The region I{r,) 


is not a fixed region. By considering the balance of the deviation 
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a’ in it, we find, by a slight modification of the argument that 


led to eq.(2.30), a contribution - [a,Gradv~ + S(a.Grad v+) 1-0 di 
J 


rt 


which represents the rate of change of "the amount of a'" contained 


in 1(P,) due: to. the tilting of the "belt" @ (T_), as =. moves and 


t t) t 
is deformed. The occurrence of the tangent vector Gradv~ in this 


CONTEXT. Ise naturals Since according, coveqn(d lls 26)% Gradv> = 


ik 


ae 
Ge 
dt 


=. 


The relation between eqs.(II.2.19-22) and eqs.(2.13), 
(2.14) and (2.24) of the present chapter will be elucidated in 


the next section. 


4) The normal moment expansion of fluxes and density integrals 


The purpose of this section is to investigate the status 
of eqs.(11.2.19-22) in the context of the present theory. Our dis- 


cussion will proceed from the exact relations 


A*(y) = | ees | a. dA 4 oat dV (4.1) 
y a X,(aV) 
o*(ay) = | JN dA + | d'°N dA (4.2) 
3V A,(V) 


that were obtained in section 2. The flux equation will be treated 


first. Let A be the arc-length parameter on L,(aV). At a point p(A) 


ic 
of this curve, the plane spanned by N and n, or equivalently by n 


and 0 scuUcs ZV) along a normal section ie Let n represent the 
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arc- length along this curve with n = 0 at p(A) and positive values 
of n in the v direction. The two parameters A,n give a parametri- 


zation r(A,n) of a region of £,(V) containing L, (av), which may 


t 
be extended to a parametrization 

R(Asn.c) = F(An) + om(Ayn) (4.3) 
of the adjacent three-dimensional space. Given (\,z), let n(A,zc) 
be that value of n such that R(A.n.C) e aV. (It is assumed that 
the geometry of 3V is simple enough -no folds- that a single n is 


so defined.) An obvious parametrization of A,(9aV) is then 


t 
BO: SG arOwe rahay (4.4) 
in terms of which the area element on A, (av) SECON =e em dare 
Wich esc) = 8p 30 . It follows immediately that 
OA 0G 
> = os eee pie => 
[Sie fa [ar Fo. fascdo.c) . (45) 
A, (dV) L,(dV)-6 


It must now be remembered that dV is a "macroscopic object". Given 
the thinness of the transition zone, a first order approximation 
of 2(r,c)N(,z) about ct = 0 will thus be very accurate. Next note 
that unless N is excessively close to Wore-n at ie the function 
n(A»t) will differ very little from 0 as zc ranges over [-6,6], 
while for any sensible choice of the dividing surface, the varia- 
tion OF an Ounee) at fixed (\,z) is slow (macroscopic). Hence, 
with great accuracy, DOGO wen may be replaced by te Reo) + 


= 


MiGAcsG 8 (J'o R)(X,0,c) for ct « [-6,6]. Hence it appears that a negli- 


on 
gible error is introduced by writing 
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on 
A, CoV) L, CoV) -6 


[2(A,0)N(A,0) +58 (ON)(A,0)] do. 
aC 


Discarding terms of order > 2 in z, this may be rewritten as 


JN dA = | fa7*NeJo + 9 (9N)(A,0) di - oN DS dh eT ae 
A, (av) L, (av) re 


> >+>>54> _ a > P 
where a = Nev, t = Nen/Ne\ and Jy and J; are as defined by eq.(2.9). 


We shall see below that this equation represents the zeroth and 
first order contributions form a systematic expansion of the left- 
hand side in terms of normal moments of J'. 

Using eq.(2.10) and remembering that N lies in the plane spanned 
by n and), we find 


on Ned) = 1 (Nefiti+ NOD) + (Jo + Jon) 


= Joe +13 
so that the first term on the right-hand side of eq.(4.4) may be 


rewritten as 


ColNGin ch =| (Ge a eh (4.7) 
L,(av) L(V) 


This represents the zeroth order contribution to the left-hand 
side of eq.(4.6). Its mathematical structure is identical to that 
of the expression used in the empirical approach to compute the 
flux (eq.(II.2.21)) so that another important link has been esta- 
blished. However in the present theory, 55 is defined as a zeroth 


order moment by i 2 Meee, and as of yet we have no equivalent 
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of eq.(I1.2.22) which defined j+ as being equal to av’. We will 
return to this point in the next chapter. 

Equation (4.7) involves actually very little information 
concerning the geometry of the surface 3V in the vicinity of Ey. As 
the finer details must be involved in an exact evaluation of the 
left-hand side of eq.(4.6), the contributions stemming from the 
higher order moments of J clearly cannot all have this simple 
form. It is interesting to find the nature of the first order cor- 
rection, which is contained in the last two terms of eq.(4.6). Un- 
fortunately their computation is surprisingly involved, the compli- 


cations arising from the term 32(4,0). We give an outline of the 


— 


derivation leading to the final result, eq.(4.15). We start from 


Q2= [dp al = dp|7}ap|2 - (a0 90)? (4.8) 
mao el San ae de 
oA 0G , OA dC L oA oC} 
To evaluate 32(i,0) we need only know 36(4,¢) and 30(A,z) to 
first order in z. From eq.(4.4), it follows that 
=. => 5 ae ee 
30(A,t) = 9R(A,0,0) + clana2R + 92R +92% 92R J 
oA 3 dT onaA dT0A ATdAATAA |n=T=0 
+ O(c?) (4.9) 
and 
90(A,c) = aA(A,0) aR(A,0,0) + 9R(A,0,0) 
AYa 3t an nia 
+ CL 82noR + (an)292R+ 20792R + 92R]| 
ag> on [dtj) an*  adconee act? |n=c=0 
Fe O(c?) : (4.10) 


The partial derivatives of R are easily evaluated. Denoting by 
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K the curvature function of uy (eq.(II.1.12)) and using the Gauss 


equation (eq.(II.1.14)), we find from eq.(4.3) 


SG OO) kf ELI) = Sa) Si 
OX on C 
VINO) = Tape ha Rees) =e aie aoa i 
\ 3 = i 5) 9U95 = 5) 9VU5 = 
ae ane eee 
92R(X,0,0) = Vev+T+S(5) mH, 92R(A,0,0) = -S(T), 92R(A,0,0) = -S(9) 
dA0N dAOT hater 


The partial derivatives of f are a bit trickier to evaluate. Let 3aV 
be described locally by the equation F(¥) = 0. Then for any 2,c 


BO GIG’ ene) = 0, from which 


HEP an(nO)eR( 1,0,0) 4 eek ,0.0) t= 0 
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fe Nen 
OOM ON) eS ea 
eye Nev 
To compute 37h(A,0), we use eq.(4.4) to get 
dC2 
Ee 
Dao Wis S) =e Cen, seme Oe ee) 
oe pa can 
Dias ae ~~ 
+ on 29 92R(A Fin peony?) 17 a7n aR(A Si Ncs Fe) 
a5} on? a2 an 
+ 9 92R(A, Oncaea) 
Ble 
from which 
POR = AC) eo) ene | (4.11) 
0C2 Gs 


An independent expression for 926(A,0) will allow us to solve for 
dG 
32n(r,0). Let us denote the covariant derivative on aV by V' and 
oc 
the shape operator of 3V by S', and write U = 96 . From the Gauss 
OC 
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equation it follov's that 


Oi Ved + UeS'(U) N 
O64 7 
he Vee! = wr ihe Ue (a2) 
dG N — 
v 
where U = Nx T (a2) 
Somchat 


Combining this last result with eq.(4.11) and taking the dot prod- 


uct with N yields 


1 Bud ee 27 NeS(¥) + t-K(v) cae N-0 82H(2,0) : 
(N-5)? er 
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AG Ame epee ined ow een eR Leence 
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Substituting all these expressions in eqs.(4.9) and (4.10), we find 


that the right-hand side of eq.(4.8) is, to first order in z, given 


On the other hand we have from eq.(II.1.10) 
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Inserting this in eq.(4.6) and developing the last term in the lat- 


ter we finally obtain 


ane Sx ah 
J'-N dA = (jaeu eerie) dh 
A, (av) L, (av) 


+ Geo ain h= wvoltinevo sedi) = (N-5)72 $'(0)-(jr+dt n)I dA 
L, (av) 
+ terms of order = 2 , (4.15) 
where 
WoseaK( i) Ky) zT a t(Neo) 2 e' (0) : 
M =Wia- 51S 1) = Ss) + Graddr 
The first term of the first order correction in eq.(4.15) (last in- 
tegral) has a structure reminiscent of the zeroth order contribu- 
tion but the analogy is lost in the other two terms. Hence, from 
the point of view of the present theory, the expression (equation 
(11.2.21)) used in the empirical approach to compute the line flux 
is a zeroth order approximation. 

The derivation of eq.(4.6) presented above was meant to 
show clearly the physical ideas that justify a first or zeroth 
order approximation. It is possible to take a more formal approach 
and obtain a systematic expansion of eq.(4.5) in terms of the nor- 


mal moments of the deviation jie The result is: 
> > eS, a 
JeNdA = J | Ad ad (4.16) 
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where a is given by eq.(2.9) and the differential operator it is 


defined by 
m 
7 zk R 
ie = ay A (0=) (4.17) 
cee Re 
with AU = ) a yn, : (4513) 


The expressions occurring in this sum are defined through the 
following equations 


agen = ee (rs = 0 for £<k) (4.19) 


=O 
NRC reece NO) 
n=0 
The zeroth and first order terms in eq.(4.16) agree with eq.(4.15). 


The fringe integral in eq.(4.1) may be analyzed along 


Similar lines. Using the same parametrization as above, let 


2 (se. or 
ran) = [aex = 
and Oc = ea OCH ee co K ea) 


It is then a simple matter to show that 


6 TCAs Gy 
Cardy. = dA dz | Gn POR Osmo la Cemse) 
i 0 
X, (dV) L,(oV) 6 
It follows from this equation that 
SED AND = [ a(a) da (4.20) 
X, (av) L, (av) 


with a = -ta; + ) CD 
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(Qe. apl4-DA , Mae Sy, é (4.21) 
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where T,H and K are the operators of multiplication by IT,H and K 


respectively and the differential operator pom) is defined by 


m 
pim) = } ] pettp_yk 
R=0u(heiD. er ae Bn 


It follows from eq.(4.21) that the fringe effect expressed 
by eq.(4.20) involves only the normal moments of a' of order > 1. 
Thus, to zeroth order, the fringe integral vanishes. 

The conclusion that may be drawn at this point is that 
when all normal moments of order = 1 can be neglected, a complete 
correspondence is achieved with the empirical approach. For not 
only the surface balance equation but also the expressions for 
A*(V) and o(aV) then have identical structures in both approaches. 
Indeed, on retaining only the zeroth order term in eq.(4.16), the 
left-hand side is approximated by the right-hand side of eq.(4.7), 
which has exactly the same structure as the expression (eq. 
(I1I1.2.21)) used in the empirical approach to compute the line 
flux. Likewise, on neglecting all moments of order = 1, the last 
term in eq.(4.1) vanishes and we are left with an expression of 
the same form as that used in the empirical approach (eq.(I1.2.19)) 
to represent the surface contribution to Dy. 

In situations where the first (or higher) order normal 
moments would play a role, the correspondence would no longer be 
complete. For, although the surface balance equations would still 
be formally identical, the expressions used to compute A”(V) and 
©*(9V) would then be different in both approaches, and the surface 


current density j appearing in the surface balance equation would 
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no longer have as simple a physical interpretation as that found 
in the empirical approach (eq.(I1.2.21)). 

To round up the connection with the empirical approach 
we shall also have to look for an interpretation of eq.(1I1.2.22) 
in the context of the general theory. We will find it when we con- 
Sider balance equations for higher order moments in the next chap- 


ter. 


5) The information contained in the normal moments. 


In this section, we shall briefly comment on the infor- 
mation contained in the normal moments. We write f(z) = a'(p+cn(p)), 
where a' is the deviation defined by eq.(2.1), and introduce the 


change of variables x = 2/6, g(x) = f(xé), so that SUP = 


Hl 
Or 
rx 
O 
V 


eee Fev” dx ie re In terms of the scalar product ( , ) 
-1 

on the Hilbert space L?[-1,1] (Reed and Simon [1972]), we have 
<g> = (g,x”). (The function g, being bounded on [-1,1], clearly 
belongs to this space.) Knowledge of all the moments Or enables 
one to compute the scalar product of g with any Legendre poly- 
nomial, and therefore determines g uniquely. Hence knowing all the 
moments soe is equivalent to knowing g itself. However the infor- 
mation contained in the full sequence (<g>) meN is mathematically 


redundant. For it is easily deduced from the MUntz theorem 


(Schwartz [1959]) that given any fixed integer mp = 0, the sequence 
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of functions a 1s*totalleinethe=spacevl"| -1.11. @1t) follows 
that mere knowledge of all moments <7, with m = my suffices to 
determine g uniquely. An immediate consequence is that if a func- 
Cicnedwer taeda lesatisties <g>, = 0 for every m= mo, then g = 0, 
and this no matter how Large mo is! This state of affairs results 
from the fact that the vectors x” are not topologically linearly 
independent in L*[-1,1], in the sense that any one of them can be 
written as a L*-converging series in terms of the others and ac- 
tually in infinitely many ways. 

The aim of the Gibbs approach is to extract gross 
features of the interfacial structure that are sufficient to des- 
cribe its influence on the bulk fields. The first few normal mo- 
ments of the deviations are adequate for that. As stated above, 
the full sequence of moments contains all the detailed informa- 
tion about the exact densities. But with mere knowledge of the 
first few moments, it is to be expected that some questions can- 
not be answered. Consider for example a flat dividing surface ore 
and suppose that we wish to compute a flux 6*(S) through a flat 


surface S panrallek to r,. As long as S lies outside the transition 


t° 


zone I this flux may be computed using the bulk current density 


be 
iy But with mere knowledge of the first few moments eee, AO 
it is strictly impossible to compute this flux when S lies in- 


side the transition zone I,- The most favorable situation is at 


the other extreme when S meets 5, at right angle everywhere, that 
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is when N = vy or equivalently Nen = 0 all along the curve SNE, 
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The more N receeds away from Vv, the less accurate becomes the 
zeroth order approximation expressed by eq.(4.7). This failure of 
the approximation is already manifest in the behavior of the fac- 
tor t which goes to infinity as N approaches n. Note that the same 
divergence is present in the empirical approach (eq.(II.2.21)). 

In that context it may be understood in purely geometrical terms. 
Consider for instance a flat dividing surface te and take for the 
region V a cone fixed in space, with axis perpendicular to Dee In 


Sy* arises from a purely con- 


the empirical approach, the term ta 
vective effect: it represents the flow of the surface excess a 
through aV, as this excess is carried along with the dividing 
Surface Dee It is clear that the rate of this convective trans- 
port becomes greater as the aperture of the cone increases. When 
the aperture is close to 180°, a small displacement of Zr, can in- 
deed produce a dramatic change of the area of the enclosed patch 


r,(V), meaning that a large amount of surface excess has crossed 


. 
the boundary 3V, in a very short time. In fact, for a given nor- 
mal velocity of Dye it 1s geometrically obvious that the rate of 
transport through the boundary via this mechanism goes to infi- 
nity as the aperture of the cone goes to 180° without reaching 
it. This divergence in the empirical approach is embodied in the 
geometrical factor t appearing in eq.(II.2.21). It is a by-product 
of the idealization implied in using excess densities concen- 


trated on a surface. In a real system this divergence cannot oc- 


cur. The fact that it is also present in the zeroth order approx- 


aval aon. (Sykes 
vale at Yo. notvadaa ade-nt \aicelgioe:  okaee 
SS? AGS THAT Stati H-eaioeoggs yv eergghdiier et eay” 
WihySethiipe) damage Chat iqen aee unt tregey et sonagye 
GT Tavita imag Vie af hod cabin st vam 2h Jeerand % =a ] ie 
at P RT G4 ie yt ee aN mite jai salt »-oamasant Vet qobte 
Sor epiustiniar rg fae eye req? nt Wentt anoo.8 ¥ sina 
ents Ghat & tioy® rears ver nia oly yRORCHS Peon Fema he 7 : 
, Seqghe. soery ps io ae th? ae. aeeeerwey 2! cds ie avid oey - 
AWTS Sus Ab ly) Onde CEP VRS ay fad sam erAT 28 M6 Agueriny 
Shad Avieneyace whi ogre ok att Aeelo-ar 3 +3 soot | 
nea J2900a ar dp tdtretvawn orl) oh aay. aamoosd Igy 


ij Aso Sey tw rs Mee Ol > “0M oF aoato 6) equrvage ant 


‘ 
ftided ‘patefoAa cet io apy et te eptats ahemes> a goubovg besb 
chasis 264 Magy a aoe ie 4 ityea pret s th ? Qh Cnor AN) 3 
-Nor. SVR B'96* ONT te Aas Dae rey ea c vwabwed unt 
bh. Stew oo) talk. wae wily: stereos a) 27 ie cr yf foofayv ‘am 
- Mae of dot meant Bed eliv (@ipweo od? Qgurias dionensie 
ip) oneg, Dupo: Drkt: oe ( PHO) 2702, « ‘oo Stas oa) 26. win | 
yal ah De raodoe el 9 MABE he foubvi gine aly int aumap ray th 2haT 2 \hY : 
ganerenoyd 8 ar A. (ASR SED) vo of Beagle |) wotaet Tes “SOM oe. 
~ genre ay! *% Neale 2K s Oni ea boop Depo? vorrentts 
=i) foanes Kutter outh bs me ANT oT 6 tel 


~eerrthe ek itiasS atlas AP Yaeiury Oate 4. 


imation of the rigorous theory in eq.(4.7) simply means that this 
approximation is no longer good enough to compute the flux when N 
becomes too close to n. The nature of the approximation is best 
understood by considering plane surfaces. The situation is il- 
lustrated in Fig.3, which shows a section of the surfaces De and 
S. The distance perpendicular to the plane of the figure is mea- 


sured by X, and the unit normal to this plane is T = nx. The 


curve Sak, is denoted by f. The position vector of a point of the 
part of S comprised between r}°and se may be written as 
ee ene ceotane y 
SE ae ie vs (549) 
and the area element on S is clearly dA = di de = didc 
SUL rane 


It follows that 


| Goetedle = | i J'(AT + ch - te¥)°N dadz 
S wee Ned 
The last term in eq.(5.1) represents a displacement parallel to 
rye Its largest norm in the region of interest is té. If the 
variation of J' along that direction and over this distance may 


be neglected, we have 


Te a 6 > =a Ae 8 
| J'*N dA =x | | J'(AT + cn) eN dA dz 
S ies Nev 
0) Sait ys = | (ja-0 4 tdq dd. 
r Nev r 
_— = _ _— 
where we have used Pais, ei eer igneand N = (Nen)n + (Nev)d. 
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Figure 3: The zeroth order current excess for flat surfaces. 
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This is identical with our zeroth order approximation. Hence we 
see that, in this case, the zeroth order approximation consists 


simply in neglecting the variation of J! parallel to 2, over a 


t 
distance té. The smaller 6, the smaller is this distance and the 
better is the approximation. The smallness of 6 ensures its va- 
lidity when N is not excessively close to n. This simple argument 
would break down for a curved dividing surface. However it follows 
from eq.(4.15) that the curvature effects leave the zeroth order 
contribution unaffected and make their first appearance in the 
first order correction. This correction takes into account the 
curvatures of both ry and S, but also keeps track, to first order, 
of the lateral variation of i: this may be seen through the ex- 
pressions of w, and 2, following eq.(4.15). It is clear from Fig.3 
that, as 6 decreases or equivalently as t increases, the lateral 
Vaniat ion OF j! may play a greater role, which means that higher 
order terms in eq.(4.16) begin to contribute, and higher order 
moments are needed. The structure of the first order correction 

in eq.(4.15) is however already quite complicated, so that higher 
order corrections are not worth considering. 

The purpose of this discussion was not to cast a doubt 
on the validity of the expressions obtained in the previous sec- 
tion. In fact, the thinness of the transition zone ensures them 
a wide range of validity. Instead we deliberately considered 


an extreme situation - probably devoid of any practical interest - 


just to show what kind of information cannot be extracted from 
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the first few moments. In conclusion we emphasize that it is only 
the approximate expressions used to compute fluxes or density 

integrals that are subject to possible limitations; the first sur- 
face balance equation (eq.(3.3)) is exact and applies without qua- 


lifications. 


6) Vector balance equations 


With a view towards the momentum balance considered in 
Chapter V, we now deduce the easy transposition of the first sur- 
face balance equation to the vector case. The idea is very simple: 
if @€ is some constant vector and A is some vector attribute, then 
@+A is a scalar attribute and as such, has densities which must 
obey the first surface balance equation (eq.(3.3)). If the sche- 
matic density is denoted by a and the exact density by a. the 


various moments a, are defined by applying eq.(2.9), a' being re- 


placed by a’ =a = a. As notediin Section 2 of Chapter II. the 
bulk current density associated to é+ a may be written as M(é) 
for some linear operator field M. The deviation between exact M, 
which we denote by M, and not M” to avoid confusion with the ad- 
HOiInts Of Maas = ML - M and its normal moments Me are defined 
in the usual fashion. Denoting by jg the surface current density 
associated to é- ae, we have from eq. (3.4) 

é 


jg = PM,(é) + SPM,(@) + @°a, Gradv’ + @*a, S(Gradv) 


Writing as usual @ = P(@), this relation may be reexpressed as 
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Jy = MP(8) + e-nT (6.1) 
where 
Wee PM)P + SPM,P + Gradv' @ a; + S(Gradv-) @ a» (6.2) 
and 
T = (PM) +SPM,)(n) + ayGrad v~ + az (Grad v~) (G2) 
with a = P(a_) and a = if et. ; 


At each point p of Dye the linear operators on the right-hand 
side of eq.(6.2) must of course be understood as restricted to 
To(z,). From eqs.(3.3) and (6.1), it follows that the balance 


p 


5 >> 
equation for é-a> 


reads 
eH ed. + DiviM?(@)+@nT] - @G 


#ineM(e) = e-ay J = 0 | 
This equation is formally identical with eq.(2.29) so that the 
steps that led from the latter to eq.(2.31) may be repeated,yield- 


ing the vector surface balance equation in the general theory: 


da. : 2Hv"a. + DivM> + Tr(SM°)n + (DivE)n 


(ok 
Gate 


- S(T) - 6, + EM (nm) - v'aD = Os ae 


This is gormakly identical with eq.(2.31). One difference is 
naturally that we now have detailed interpretations for aes G. 


and M> and Tf (eqs.(6.2) and (6.3)). 
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CHAPTER IV 


BALANCE EQUATIONS OF HIGHER NORMAL MOMENTS 


In Chapter III, we obtained the first surface balance 
equation which was formally identical with that used in the empi- 
rical approach. This provided us with a detailed interpretation 
of the surface fields used in the empirical approach in terms of 
normal moments of the deviations. It is interesting to know that 
the higher order moments also obey exact balance equations of a 
Similar nature, though not involving the bulk fields. It is the 
object of the present chapter to derive these equations and dis- 
cuss the terms that enter them. Despite the fact that these equa- 
tions, unlike the first surface balance equation, have no equi- 
valent in the empirical approach, the first of them will provide 
us with an interpretation of the transversality condition (eq. 
(II1.2.22)) in the context of the general theory. This will be 


discussed in Section 3. 


1) Balance equations for individual normal moments 


The goal of this section is to obtain balance equations 
for the normal moments a2 m 2 a, defined by eq. (11l.2,9) 
A cartesian frame (ixed in Space is introduced, whose 


origin coincides with a point of the dividing surface aa at time 
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to and whose z axis is parallel to the unit normal no at that point. 


Rromseqss(iil.2.10)eand (111.2.8)weshave. at time tp. with m = 1, 


> ee 
ie (0) 5 figrd _, (0) = | Zs (ztieeMandz 


The limits of integration have been pushed to infinity in view of 
later simplifications in the notations. The jump in J} at 6 
makes no contribution in the integration by parts, thanks to the 


factor 2” (m>1) . Finally the fact that ad3 is not defined at 
OZ 
0 is obviously irrelevant. With the help of eq.(III.2.3),the last 


equation may be rewritten as 


ae ga adi , adi 
ee en) me Ox Oy °n (1.1) 
We first concentrate on <2 > To bring out the moment a> we 


write 


eo Ss = | Z eae (zine) dz = d_| z"a'(zMo,t) dz | s 
t= 


Gin) 
Here again thanks to the factor z”, the jump in a' makes no con- 
tribution. At time t the dividing surface cuts the z axis at 
Ze=eLc) cain pati Culdt eZ ty) = sOumniten Q(t) defined by Q(t) = 
ae . c a Ale Me 3 
Z(t) ig » 1b 1s easy to see that Q(t,) = Vv ny, where Vv) is the 


normal velocity of ry at cr The difference no - n(Q(t)) will be 
0 


denoted by An(t). After making the change of variable z = Z(t)+z 
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in the last integral, we get 


m 


[2a (2R,t) dz) m: is [ota c(zee let] dz 
N=O(M-A) tA! 


tee) ='©O) 


which on account of Z(t,) = 0 and Z(to) = v", gives 


= mvta =) + d_| 2 at iLee) ns allel dz 
t=ty uw t 


tata 


=0O 


(1.3) 


Now at any point q not on %, we may write 


t 
a'(q+th,t) = a'(q,t) + heva'(G,t) + [h[2R(q,tsh) 
where the function R is bounded on any bounded domain of variation 


of its arguments. Applying this to the integrand in eq.(1.3) gives 


| ee iec\tigst) do = | ea 1 Zngecneean, cl de 


—0O -CoO 


= | (a (inert t )tcan-Val( 2hoten,t)2l ean Ri cng renst scan) | dz 


m+2 


a (Q(t) + An(t)«[<va enon 3 + | An(t) Pls R(Zno+tn,t3cdAn)dz. 


=0o 


4) 


> > 
BUcetromscae( fl.) 9), 1) fol lOWSathatedsaneUutyye= da (0,to) 
dt” dt 
while, from eq-(II.1.26),d An(t,) = Gradv* (6), and finally the 
dt 
time derivative of the last term in eq.(1.4) vanishes at to since 


= => 


on account of |An|? = 2[1 - norn(Q)J, both |An|? and its time de- 


rivative vanish at ty. There results the following expression for 


a di dL A t 
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The last term in eq.(1.1) will now be rewritten in a mani- 
festly invariant form. From the definition of the directional deri- 


vative, it follows that for any U « Ts(Z,) 


D-<J'>_(p) - d_ | MSE le) + om(Ece))1 a : 


de). =e 
>. x =e at Ss = 
where C is any curve in Z, such that C(O) == prandecG. 0) = sureaine 
e 
right-hand side may be rewritten as 
“om 15 (2.2/2 e ' 71 
| C Lhe aenye 1 (Beon(p)) dz = <D.J'> - seas hore 


SOuthac 


ae 
<D5d ee D-<J a <Dg(g) J" met (126) 


Introducing the covariant derivative Tak defined by eq. (11.1.4) 
and using eqs.(III.2.10) and (II.1.14), the first term on the right 
may be rewritten as 
' _ rod - | _ a = oul > 
D-<J a Vide O10) aU (S(j,,)+Grad J) J (anlar) 
Introducing likewise the covariant derivative se which, for 


each fixed q not on X,, iS a linear operator on IR? defined by 


£ 
4 > >> . A 
(DJ') eh) : DeJ'(q)> the second term in eq.(1.6) may be rewritten 
as 
> os ms 
<Dg(q) 4 > cna <(DJ JS ; 


It. follows that 


>’ 


Ou; dU5 _ ihey ek — - 1 Tat 1 
coil 4 S?> = Divi, - 2Hi. + Tr<P(DJ')S> 


OX dy m+1 


where use has been made of eq.(II.1.13) and where P is the pro- 
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jector on Tay) 
Combining this equation with eqs.(1.1) and (1.5) finally 


yields the requisite balance equation 


2 we {Os Hien eal ih . 
an 2Hj. PD AV ea Chace, ma av “Alec + Gradv -<Va ave 
+ Tr<P(DJ')S> = 0 (Ves) 


m+1 
which,needless to say, applies only when m = 1. Some comments 
about the last two terms are in order. Clearly, the fact that Va' 
and DJ' are not defined on the dividing surface is irrelevant in 
computing the moments. (Derivatives are not meant in the sense of 
distribution theory.) Next we shall consider the size of these 
terms. The first one may be rewritten as 


dL ' = 1. 1 
Gradv-<Va aes Grad v+P(<Va an 


= i, 
= Gradv-<P(Va ) erect 
Whence this scalar product actually involves only the tangential 
part of oC aa . Because variations of the deviations along sur- 


faces parallel to x, are comparatively slow, it follows that 


c 
Grad v"+<Va'> should be small compared to the other terms 


td 

of order m-1 and m appearing in eq.(1.8). This would clearly not 
have been the case, had the normal component of dian, been 
involved. Likewise, the term Tr<P(Dd")S> ought to be small 
compared to terms of order m-1 and m, thanks to the factor 
S which insures that only tangential derivatives of J! are 
picked up. 

An apparent drawback of eq.(1.8) is that the last two 
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erates a series converging to P(<Va'>,): 


P(<va'> eSaGaada (lad) 


b) Reo R+A 
But under usual physical conditions this series will converge so 
rapidly that the first term will already be an excellent approxi- 
mation. In any case, eq.(1.10) provides an explicit expression 

of the remainder at any order of truncation. 


A similar recurrence relation can be derived for 


<P(DJ")S>, . From eqs.(1.6) and (1.7), it follows that for any 
Ure Ts(Z4), 
<P(DJ') P>, (u) = P(<D-J'>,) 
Ieee al a = 
= Ft JpS) (a) + <P(DJ )S>, (4) 


from which 


<P(DJ')P>, = V5, - ins + <P(DJ"\S> Gn) 


Multiplying on the right side by S gives the requisite recurrence 
relation 
41 AG 74 
<P(DJ )S>), = (Vi, J,S)S + <P(DJ )S>) 48 : CRESS) 


Upon writing 2, =p - aS. it follows that for any integer 


AAW a (1.14) 


The comments that were made about the next to last term in 

a > 
eq.(1.8) apply verbatim to the last term Tr<P(DJ | Soe . Thus 
form eq.(1.14) we have, with a high degree of precision 


Tr<P(DJ')S>, ~ Tr(2,S) = Tr(SVj,) fe FS) 
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Trt (S+ 2HI)V5,1 - 57, (4H? - 2K) 
: TLSV J, ] + 2HDivj, + 5,(2K-4H?) 
Div S(j,) + 2HDiv J, + jj,(2K-4H? ) 


where in the last step we have used eq.(II.1.9) together with the 
identity Divs = 0, which is equivalent to 

Diesen =u Grad He. C115) 
The proof of this equation will be given in the next section. 
Better approximations could be obtained if necessary from eq.(1.14). 
The procedure for obtaining balance equations for normal moments 
of any order is now clear: one starts with the basic relation ex- 
pressed by eq.(1.8) and uses eqs.(1.10) and (1.14) to reexpress 
the last two terms in eq.(1.8) if necessary. One thereby obtains 
a hierarchy whose truncation can be based on clear physical argu- 


ments. 


2) Balance equations for combined normal moments 


The balance equations derived in the previous section 
(eq.(1.8)) involved terms (the last two) which were not expressed 
directly in terms of the moments used up to that point (a), jp and 
ip ). We showed how this could be remedied at the expense of very 
mild approximations. In this section we shall see that the balance 
equations for the normal moments of order = 1 can be combined 
three at a time to obtain exact equations very similar to the first 


surface balance equation of Chapter III, and involving only moments 
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of the type ap, jp and ie 


From eq. (1.11) together with -2HI = -S+S and SS = SS = 
=Kl (see eqs. (111.2720) and (111-2.22)). we have 


PeVal> = 2H-val> + K<Va'> 


k+1 bro) 


- ) s*(Grada,_.) - y S***(Grada 


R+A 


~N 
+ ) SS (Grada 


R+A+1 


= Grada, + S(Grada, ) ‘ (28 1) 


Likewise with 


M, = <P(DJ')S>), - 2H<P(Dd')S>) + K<P(Dd')S>, (2.2) 


we have, using eq.(1.14), 


co 


Gage 
-)2 
nel Raat 


A+1 


, fio) ee iw OS 
Meee) CS S +) Q Sees =) 
k ee R+A = R+A+1 jae R+A+2 


= QS: + % ,495 = 2,5 - KQ (23) 


k 
We calculated Tr(2,S) at the end of the previous section, and 
ee aa ee all. = 3. So . orl ieee 
Tr(2, 4) = AVA tes 4 de>) = USE AW . Combining these 
results we find that 
lie = Div£S(j,)1 + 2HDiv 3, = KDI 

ON oak oullt 
TuCCieee eel 1 + ANI 5 peed) 

Let us define the following combined moments in terms of the 


normal moments defined by eqs.(III.2.9,10) : 
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When m = 0 , these are the surface excess densities (see eqs. 
(III.2.11,12)). We now pick three of the balance equations (eq.(1.8)) 
of order m,mt+t1,m+2 (m > 1), multiply the second one by -2H and the 
third one by K and add up all three to get (using eqs.(2.1) and 
62525) 


Is -LS ale att 0. . = 
da. - ens ~ Za jdt - a aod kK + Div J, ~ cHDiv 4 
c e dt 
5. oS S fh lll: dk ot 
+ KDI Jee =o ma _4v = ee) - 2H(m+1) (a, v - Jn) 
+ K(m+2) (a, ae red + Gradv'+(Grada_ + S(Grada__)] + DIM ie =) 
(256) 


Now we shall show at the end of this section that 


d-H = (2H2-K)v~ + #Lapv- (Wapwe= BaviGrady (a2) 


and 


dk] 2Hky = Div eS(Cradv-) 1) (2.8) 


Using the self-adjointness of S we have 
ii = Se ae ik 
Grad v*[Grad CW ee: S(Grada__ yale Div(a_ Grad Vex) Qe LAP V 


; Ps Jt : Pe iN: 
+ DivCa_,., S(Gradv ee ano Div [S(Grad v )] 


Inserting this together with eqs.(2.4), (2.7) and (2.8) into 


eq.(2.6), the latter reduces to 
LS LS ba ee @ iz L ps i 
aa, - 2Hv at Diab ra SIR cr) + ana Grad v ‘ an.oo (Grad v fa 
Sl ees 
=o 2 mas - Jae) 0 (2.9) 
(m > 1) 


This equation is exact and has a remarkable similarity of struc- 
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ture with the first surface balance equation, eq.(III.3.3). This 
Suggests that the same technique could be used to derive the first 
surface balance equation. Indeed, eq.(1.8) may be shown to hold 
also for m = 0 provided that the term ma vo . de be replaced 
by - fied ie al ved Enron - av » and the same calculations as 
above may be repeated, leading eventually to eq.(III.3.3). But 
this procedure would only establish the first surface balance equa- 
tion as a necessary condition for the formal balance equation 
(eq.(III.3.1)) to hold, whereas the derivation given in Chapter III 
shows clearly that it is both a necessary and Augpictent condition. 

We will now give the proofs of eqs.(2.7) and (2.8) and 
next of eq.(1.15), which were used in deriving eq. (2.9). We take 
ae at a fixed time t as a reference configuration. If U;,U. are 
two tangent vector fields, linearly independent on each tangent 
plane, we have 

2H = Uy > S(u,) + Up S(U2) (2.10) 
and K = [Ur + S(t.) ICU - S(G2)] - CU + S(Ui)ILU2 +» S(t2)I (2.11) 


where (Us su.) 4s the basis reciprocal to (U,,U,), that 1s 


u., Up 6 g? so that 
Uy; = Uoxn U> = nx, (2212) 
Q Q 


where Q=n+(U;xU»). Now we let U, and Us be "normally convected" 


on 2a according to the defina tionsof secticn ItoteChapter Tis 1 


is then a simple matter to show that 
d'S(v_) = Ve Gradv~ + S(i,):Grad vn 


dt c=0 Qa 
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while from eqs (II.1.24) and (II.1.25) 


a = -vS(U_) + U.+Gradv n 
dt T=0 2 ” 
and 
d= = - 2Hv~ 
dt T=0 


Combining these results with the normal time derivatives of eqs. 
(2.10,11) and doing some algebra, we get eqs.(2.7) and (2.8). The 
calculation is much simplified if U,, U2 are chosen to be eigen- 
vectors of S on the reference configuration ae which is always 
possible. 

We now come to eq.(1.15). Although it looks very simple, 
its proof is not so straightforward. The formula is likely to be 
known to some differential geometers, but we could not trace it 
in any of the several standard references we looked through. The 
crucial ingredient in the proof is the Codazzi-Mainardi identity. 
Let p be a point of a surface © and U;, U2 be two orthogonal unit 
tangent vector fields defined on a neighborhood of p in = with 


U2 = nx, . From eq.(II.1.8), we have 


a Divise=seuq ((V, S) (ap) s (VeeS) (us) ie: 231 3)) 


1 U2 
The first term on the right may be developed as follows: 


Uy ° (V5 S) (an) = UieV CS(ui) = Ui *S(V> ui) 
= Vy Cui S(t )I S 2S a av. Ue 
SU, Grad 2H Vella -S( Us) 1) =e 2S\U,) <Vonus 
Ui Uy 
=U, °Grad 2H = S(U2)+Ve Uoees Ue Vous (Ua) sl 
1 Uy 


= 2S(ti)*V5 ta 
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= U,° Grad 2H - U2? V7 (S(t) J 

oa S(Up) +( 

(2.14) 

In the last step, use has been made of the Codazzi-Mainardi equa- 
tion 


Vy S(¥) - Vy S(X) = S([X,Y]) 


where [X,Y] is the Lie bracket of the tangent vector fields X and Y 
(Hicks [1965],p.8,28).The second term in eq.(2.13) may be rewritten 


Wir (V5 S)(2) = Uae {V5 CS(G2)] - S(T, Ue)} 


where in the last step, use has been made of the identity 

V; ue = Vq,u = fU;.Uoh (Hicks £19654, p-26). Inserting this. 
together with eq.(2.14) in eq.(2.13) yields 

Ui DivS = uy*Grad 2H - 20S (U2) *V5 Ue + S(ti)-V5 tJ : C27) 
We will now show that the expression inside the square brackets 
vanishes. Since U2 = hx Ui, we have 


V, U2 = PLD; (nx Gj) = PL(D 
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and S(U2)+V5 Us + S(Ui)+V, Uy = [-n x S(Hi2) + S(t) 1+V, Uy 


But since U,; is everywhere a unit vector, it follows that Vy, 4 
is perpendicular to U,, that is V5,4 = aU, and the last expres- 
Si0n may be rewritten as 

[-nxS(U2) + S(y) Jeaila = af (NxU2)+S(G2) + $(U1)+O2) 


of-UieS(U2)+S(ti)*t2] = 0 


from the self-adjointness of S. Equation (2.15) thus reduces to 
UieDivS = U,*Grad 2H 
The same relation obviously holds for U2, so that 


DIV Ss = Grad 2H. 


3) Application to the transversality condition of the BAM theory 


In the previous chapter, we have provided the bridge 
between most of the equations underlying the empirical approach 
(eqs.(11.2.19-24)) and the general theory based on deviations. 
Equation (II.2.24) was given a rigorous exact analogue whereas 
eqs.(I1.2.19-21) (especially the last one) were recovered in the 
context of a zeroth order approximation. However one last gap 
remained to be filled, eq.(I1.2.22) 

je eaves (Sis) 
What we found in Section 4 of Chapter III was that the zeroth 


order contribution in eq.(III.4.16) gave 
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This has the same structure as the line integral in eq.(II.2.21) 
except that in the general theory i; is, according to eqs. 
GlIite2. 9310) seciven by 

i (ee oe rie (3.3) 
and not by an analogue of eq.(3.1). Since eq.(3.2) was obtained 
in the context of the zeroth order approximation, we might expect 
that eqs.(3.1) and (3.3) could be reconciled within the same con- 
text. This is confirmed by the balance equations derived in Sec- 
tion 1. Indeed the balance equation for the first order normal mo- 
ment obtained by setting m = 1 in eq.(1.8) reads 


dee eet ice Di Vie moe x anv 


oO 
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puGtadive-va er Tr<P(Dd')S>. = OC, 


In the zeroth order approximation we neglect all moments of order 


IV 


1 . The above equation then reduces to 

ie See (3.4) 
Since in the zeroth order approximation the surface excess a> is 
set equal to ao (see eq.(III.2.11)), we have reconciled eqs.(3.1) 
and (3.3). The whole set of equations (II.2.19-24) which expresses 
the structure of the empirical approach has thus been given a firm 
basis in the context of the general theory. Equation (3.4) agrees 
with the transversality condition derived in the theory of Bedeaux, 


Albano and Mazur, as may be found for example in (Vodak [1978], 


eq.(13)). 
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Having obtained exact balance equations in terms of nor- 
mal moments, our main concern has been to interpret the empirical 
formalism in terms of normal moments of the deviations. This could 
be achieved in the zeroth order approximation. It is to be ex- 
pected that higher order approximations could be used to handle 
Situations where very high curvatures are present or where the 
interface is not so thin. Buff [1956] has thus studied curvature 
effects in equilibrium . The theory developed in this chapter 
and in Chapter III could provide the basis for such an investi- 
gation in non-equilibrium situations. However at this early stage, 
we restrict our attention to features present in the empirical 


framework described in Chapter II. 
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CHAPTER V 


MOMENTUM AND ENERGY 


1) Mass and momentum balance 


Let us consider an n-component system without chemical 


reactions, with schematic (i.e. bulk) and surface excess mass den- 


Sities eae. and Pepa We recall from Chapter III that 
on = po - 2Hor + Koo Cit) 
: Ch Qy 4 
with Pn = <(on) ek 


: see 4 =n all 
The associated current densities are denoted by J,...,J andj ,... 


=n 
j where 


Ge jn ee Sji) + Grad V4 o55( Grado) (ls) 


and is = p(s”) 2 Aas 


The total densities are 


n 
Cao) 
Ci=a 


sae Gin) 
all 


i715 


n 
) oe peel G = 


se palais 

a 
Py > Ge) ds Ge = 
a=1 al a 


a 
The general balance equation of Chapter III together with the usual 


bulk balance equation become 
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as 


op - ¥V°Gu=70 (elistey. 
at 
d*p. - 2Hv+p, + Divg+ [G-n- pv'D =0. Cha 
dt 


Equation (1.5) lends itself to a simple intuitive interpretation. 
On rewriting the last term as [o” (w’*n-v-)7, we see that 
Cane ov I may be interpreted as the rate of transfer of 
constituant a into the interface from the + phase minus the rate 
of transfer of constituant a from the interface to the - phase. 
This,together with the lateral transfer within the interface em- 
bodied in Divj- » gives the net rate of change of local adsorp- 
tion of constituant a. In the zeroth order approximation, only the 
first terms on the right-hand side of eqs.(1.1) and (1.2) are re- 
tained. 

The bulk total mass current density é and the bulk mo- 
mentum density are one and the same. If II denotes the momentum 
flux density tensor (Landau and Lifshitz [1966]) and F the body 
force per unit volume, the bulk momentum balance equation reads 


> a 


i se os (1.8) 


ot 
As before we ignore the subtleties associated with diffusion: par- 
tial stress tensors,etc.(Kehlen and Baranowski [1977]).The momen- 


> 
tum surface excess density G. 1s 


G. Pa Ca neh Gua Kes (1.9) 
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Note that its tangential part P(G.) is not equal to the vector g 
ofimeq. (1.3). sFor onedefining 
ee =e 

i, = Ln = P(<G'> ) (1.10) 
we have from eqs.(1.2) and (1.3) 

G@= dor Sta) + oiGradve 4 osSiGrady)) (eddy) 
whereas, according to eq.(1.9), 

P(G_) = Go - 2HGi + KGa . (iea2) 
This is in marked contrast with what we had in the empirical ap- 
proach (eq.(1I1.3.9)). The difference disappears again in the zeroth 
order approximation for the right-hand sides of eqs.(1.11) and (1.12) 
then both reduce to go. 


th 


Let us denote by le the m~ order normal moment of the 


deviation II' of the momentum flux density tensor: 


ley] <i> 
m m 
and write 
> 
g- = Noels 
m m 


Then it follows from eqs.(III.6.1-3) that the momentum surface ba- 


lance equation is 


d*G : 2Hv"G. + Divi? + Tr(Sm>) n+ (DivT, ) n 


dt 
-S(F,) - f, + cm(®)- v'Gn = 6 (1.13) 
where 
> = Pj) SPlyP 2 Grady eld, + o(Gradv)ergs sel 14) 
and 


ry 
i 


1 (Pies Spiny ich eregaGradve + goS(Gradv-) . Gi 5) 


LZ 


hee 


ie ws 
ids °F 
: b te pa" of ures 4 

est) t) | a 1 ONS > it oe 
=tis faatytives 5 anh a? jhe Sy iit et sin a 
dyeawkoort nt Hoge thopaaaelh Sash aeT (1 B.8, mad 


{i} : ere gs 


(ett ygine ial), ene Vol ithe neigh ta “ae not p sat eerrgne 
7 | 
) -of of soubor pete 
dhe VA Senn Teron “148% ond ai Ye atonal 2p dad 7 


eo 
amet et fart vod? inane att ta '" wotde ¥ ap 


~ iv = ae 


«(ean S 4, 

B, er ae 

“ne anabina mati sit hee ee last Ply #8 ov? awok lat vi" 
a optieaps.s 


i Cah ea ; ered ip  *paith« pPvns. ie: 7 
T = - + 
(ei, 1) A tev lat rh. 4 - (he 


te.1) °° sed wheat « Pe he a 4m. wats = : 


oo, 
' ~* 


Equation (1.13) is formally identical with the abstract empirical 
momentum balance equation (II.3.11) provided that g+p vin is writ- 
ten as G@ in the latter. But the analogy is not complete. First in 
the empirical approach, we had @> = G+p-v-n (see eq.(II.3.9)) so 
that P(G?) was equal to g, or in other words the tangential com- 
ponent of the surface momentum density was equal to surface mass 
current density. Intuition based on experience with bulk densities 
would certainly suggest such an equality. However on comparing eqs. 
(1.11) and (1.12), we come to the conclusion that it does not hold 
in the present context. Not Surprisingly intuition is vindicated 
in the zeroth order approximation because P(G°) and g are then both 
equal to go. 

Secondly,it follows from the exact balance equation 
(IV.2.9) that in general Gon is not equal to ae as was the case 
in the empirical theory (eq.(II.3.9)). Yet another difference is 
that the relation T 


LG 
the basis of angular momentum balance does not hold in the general 


= v'g obtained in the empirical approach on 


theory (compare eqs.(1.11) and (1.15)). 
Now it will be recalled that the relations n-G = ov" 
and Ty = vg were essential in going from eq.(II.3.1) to eq. 
(11.3.31). Consequently the latter intuitively appealing form of mo- 
mentum balance equation has no exact analogue in the general theory. 
It will come as no surprise now that complete analogy is 


however restored in the zeroth order approximation. In this approxi- 


mation, all normal moments of order = 1 are neglected. The surface 


118 


aban oa 


7 erty ee 


~wi09 cardindiiest wii nee 
ashe Suits bk Taupe aon \ a | 
pevetensi tind igtw asiitceanorg ene 
ape pnt reanca: lo. TayaWi 9] eel gh ie ieee 
Bilin gor ead di sare qiser ro ai wa anny oe (801) 
hakeorariv, 21 owisiuint etait anigd fom axainos 
ot nary 4x6 2 ian (78)9 chvatod abasainoreas sh diores on 
rT) at 
noPeehe siacien toes ett dot? geet tot Ih, ylemoose 
rand.282 Siw. v ie od. 1EORe Sow 2 ite a fwianep wt Jedd te.$ vi) 
af eyhare hha obazdin, 28Y ee. e4P ded enoeet (iors tqms oft ab 
an (bode tynbd fine gtd ie bantetadep we a pottaten att ana 
tarovey addon ole gon var eon, c out norvens 9 arene. 
| (oat. See TY, ape seine) 

4% = eH eHtohlaty Gitdet?. bafroagy ad bite at wow | = 
19 od (C.-T atey of Orrangebe ate g'¥ = = shone 
“ah Lo. NO? cir doges, /favedaienl sedate! aay yTimaypaened (1 SeTY 
aonny Ferorndy sti pt sipel bre tone wo ant not rnpe sant 6d, ) 
di ete Bi shetums Tate wht oeierws on 26. oineg. Pit a. 
oberg aft a NGM NiNUA gS. Went codes o@h We beret eet 


esate sat atiatgan aM i Ae-agieen Fs smmeonhianise fe 


7) 
er 


rv _ Di 


oie) 


excess mass density e. is thus replaced by p% while j” is replaced 
by 39 (see eqs. (1.1) and (1.2)). Likewise eqs.(1.9),(1.14) and 


(1.15) are replaced by 


ae Ga (Gini) 
f= lye (197) 
Ep PIlg(n) . (1.18) 


In the zeroth order approximation, eq.(IV.1.8) also 

provides us with the general BAM transversality condition 

jo = agv" (1.19) 
which we will now put to work. Since G is the current density 
associated to p, this condition gives us 

NeGo SS avews (1220) 
Combining this with eqs.(1.16) and (1.10) we find that in the zeroth 
order approximation 
G. Che] Gone Oavan (ion) 
which is in agreement with eq.(II.3.9). We will now see that from 
the transversality condition and from eq.(1.18) we recover the re- 
lation Ty = vg which was obtained on the basis of angular momen- 
tum balance in the empirical approach. Let @ be a constant vector 
and consider the scalar density é+G with associated current density 
m(é). The corresponding surface density, in the zeroth order approx- 
IMac LOM sn ais Bb while the zeroth order normal moment associated 
to the current density M(e) is Ip(é) = <II'>)(@). The general trans- 


versality condition (eq.(1.19)), for these particular quantities, 


takes the form 
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Nellg(@) = Gov. (1.22) 
Now we shall assume that II and I, are symmetric, so that so is Il), 
from which nello(@) = @*Io(n). Since e is arbitrary, it follows that 

Mo(n) = v-Go (1.23) 
which, combined with eqs.(1.18) and (1.21) gives 
ee vigem® (1.24) 
The balance of angular momentum is used here implicitly through the 
assumed symmetry of II,, which enables us to go from eq.(1.22) to 
eq.(1.23). Substituting the zeroth order expressions for G. and ry 
into eq.(1.13) gives the momentum balance equation in the zeroth 


order approximation : 


d'(Go + Oov-N) - 2Hv(Go+ povn) + Divi? + Tr(Sm)n 


+ Div(v'g,)n- S(v'g,) - fo + O(n) - viGD - 0. 


We may now of course rewrite it as 


poDV = DivT> + Tr(ST°)n+ fh + L(T- ow, ew.) (n)D (1.25) 
Dt acura 

wnere 
V -yoe Vv he. Cov =e 
T =-0 +povev . (1.26) 
We =w - Von Lt 

and D_ is defined as previously in terms of V . 

Dt 


The zeroth order expression for the surface stress tensor 


T° follows from eqs.(1.26) and (1.17) 
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= - PL<(-T + pw@ W)'>oP + poVe@V 
or 
T° = PToP - <(pW,@Wy)'>o + PVE V (1.27) 


where ¥ = P(w). The absence of any curvature term in this expres- 


sion is a consequence of the zeroth order approximation (contrast 


with eq. (1.14)). Naturally at each point p <« 2,, T° is a linear 


iE” 


operator on Tale ), so that the terms on the right-hand side of 


t 
eqn Gla2/)e.(oreforithatmattersot eqse@iadd pile 7 cand @iaco)) 


actually stand for the restriction to Tale of the corresponding 


,) 
linear operators. Our final expression for the surface stress ten- 
sor shows that it arises from excess forces (PT,)P) and from a rather 
peculiar kinetic term, to which we shall return when we take up the 


question of kinetic energy from where we had left it in Section 4 of 


Chapter II. 


2) Relative velocity residuals 


In this section, we return to the question of energy with 
which we had closed Chapter II. This quantity, as well as related 
ones which will be considered below, introduces some complications 
because of its nonlinearity in other quantities. For the first time 
we meet with expressions that, even in the zeroth order approxima- 
tion, may be at variance with those used in the empirical approach. 
Interestingly enough,the quantities or relations involved are precise- 
ly those which we could not justify on the basis of general prin- 


ciples in our attempted deductive approach to the empirical for- 
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malism in Chapter II. Interesting too, is the fact that all the 
potential discrepancies can be reduced to terms having the same 
basic structure and intimately connected with the behavior of the 
velocity distribution in the transition zone. In view of the com- 
plexities involved, we restrict the discussion to the zeroth order 
approximation at the outset. 
We start with the kinetic energy excess. The tangential 

Surface velocity field is defined as before by go = OoV where 
Ge Pica le Moreover, in the zeroth order approximation, the BAM 
transversality condition holds so that we may write NG = Dov. 
We then define a non-tangential velocity field on os by the rela- 
tion 

Weaver yon: (Pay 
and extend it to the neighboring three-dimensional space by de- 
fining it as constant along the normal at each point of the divid- 
ing surface. Finally we define a relative velocity Wi, in this re- 
gion by 

w,(b+on) = w(p+en) - WB), Pe 2, 
or briefly 


iy = Geis (a) 


The velocity w is the bulk or "schematic" velocity. In terms of the 


exact velocity w , we define another relative velocity 
weewe-V (2.3) 


and denote the deviation w. - Ww, by Wy We are now ready to start. 
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Taking the zeroth order normal moment of both sides and denoting 
by Ky the zeroth order energy surface excess AOE. we get 

Ko = 400V° + 4<(pwi)'>o . (ey) 
Thus, on top of the expected 4o,V*, we find a residual term which 
can be either positive or negative. We note incidentally that the 
last two terms in the zeroth order expression for the surface stress 


tensor can be combined into an expression similar to the above re- 


Sidual. Using the same trick we write 


(ow,ew,)' = o*(w'-V)e@(w-V) - o(w-V)e (w- V) 


+k sS ss kK>* 7 a Se aR *>* 
= O0W @W -P9WEW-0 WEVt OWeV - VEO Ww 


ok 
+ Ve ow + (0 


from which 
= > Afi ke =e > 
<(ow@W)'>o = poVeV + <(pw,2 w,)'>o ae 
Combining this with eq.(1.27) gives 
aia? =] a2 ono Pk (2.5) 
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Ignoring external fields for simplicity, the balance 


equations for the exact kinetic and internal energy densities are 


known to be 


3 (p"w?) + VeCow?w - T (W')1 + Tr(T vw) = 0 
Ate 22 Z 
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Tie TeCu we eel - Tr(T vw.) = 0 
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where, as we recall, the notation T, for exact T is used instead 
of T~ to avoid confusion with the adjoint of T. In the zeroth order 
approximation, the surface source term for internal energy is there- 
fore given by 

<Trl (TVW) 'I>o = TrE<(TVW)'>o] 
Let us see how this compares to the source term considered in Sec- 
tion 4 of Chapter II. Still using the same trick, we write 


(TVW, )' = Tue a0) 6 eis 


(Tyw)' - T'W 
from which 
> A => 

Tee <Ch ww) oel ="1 (<7 Wee) = TES) 0d : (2.6) 
In order to write the first term on the right-hand side in terms of 
known normal moments, we must find an explicit expression for Ww. 
To this end we introduce a parametrization r(&!,&2) of Zr, and ex- 
tend it to a parametrization 

piel e? Fe fest ed eS 

R(E Ae i) = UE ANS ) teGN) 
of the neighboring three-dimensional space. The extension of V to 
that space is of course given by 

oe > 

ieee ga) eRe a) = Wipes a2) 

a a nt ee os ; : ; 
so that (nev) V = 0. Now restricting the discussion to one partic- 
ular point p of Dy» we may choose the parametrization r(&!,=2) in 


such a way that, at a point ie ar is a unit eigenvector o, of the 
age 
shape operator Ss. It then follows from eq.(II.1.10) that at 


point Dp 
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where KY 1s the principal curvature in direction ey: Now we have 


on the one hand 
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Comparing these two relations, we conclude that 
(WV) (1- cS) = DV, 
from which (since (nev)V = (VV)(R) = 0) 
vv = (VW)P = (DVo)(I-cS)"'P. 
BUtEUSIhd eds... 1iel210) sand (ill. t4)eett followss that if U 1S any, 


vector in Ta(24) 


so that finally 
ee vive veSMmein rus vecrady- 1) | teres aay 
where, as is easily checked, 
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with A(t) = 1- 2HZ + Kz* . The z dependence of ww comes through the 
factor (I1-zS)~!. To zeroth order in c, we then obtain 


Tr(<T'VV>0) = TrlTo(2+ neU)PI (2.7) 
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where 


and 
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Now we have 


Tr (heQe m= sir Pier) 


= Tr(T°2) + Tr(PWPQ) (2) 
where 
W = <(pwW,@W,)'>o (2.9) 


and use has been made of eq.(2.5). From the zeroth order relation 


IIo (n) = We Go 
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It follows that 


Tr(T, neU P) = Tr£€(T (nN) @ U)P] 
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[S(V) + Grad v-J*W(n) 
Combining this with eqs.(2.6-8) finally gives 
TrL<(TVW)'>o] = TRLT?(Vv¥ - v'S)] + Trt PWP(VV¥ - v'S)] 


+ [8(V¥) + Grad v7]-W(n) + Tri<(TWW,)'>01 (ORI 
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as the surface source term for internal energy and -1 times that 
for the kinetic energy. 

The other surface densities and current densities related 
to energy may be treated in similar fashion. We only quote the re- 
sults. The convective kinetic energy current density is 

seer = tooV2¥ + 3(TrW)V 


+ PW(V) + <(pweii,) '>o Cae 
2 


and the corresponding conductive current density is 

- PL<(T(W))'>o] = - T9(¥) - PLW(V) + <(TCW,))'>0] . (2.13) 
Finally the internal energy convective current density is 

PL<(uw)'>o] = uov + PL<(uw,)'>o] . (2.14) 

On comparing eqs.(2.3) and (2.11-14) with eqs. (11.4.3) 

and (11.4.9,10), we find that the general zeroth order theory em- 
bodies additional contributions, all of a similar nature, to the 
energy-related surface densities and current densities. It might 
be thought that the presence of these residual terms is a pathology 
peculiar to a Gibbs type of definition of the surface excess 
densities, and that some other kind of averaging would wipe them 


out. This is however most unlikely as their source is a nonlinear- 


ity that cannot be eluded. The consistency requirement (eq.(II.4.5)), 


which, to our knowledge, has hitherto been overlooked, already 
points to a difficulty even at the purely empirical level. This 
problem will require further investigation but, at the present 


time, we can only leave the matter as it stands. 
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